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Energy Losses of Small Rotors * 


By Donald J. Kenney, Ph.D. 
Research Institute 
University of Detroit 


ABSTRACT 


A free body, rotating in empty space, is predicted to lose a 
small amount of its energy to the rest of the universe. Because of 
the abstruse mathematics, founded on uncertain cosmological as- 
sumptions (flat space, curved space, anisotropic space, etc.) which 
predict this energy loss, an experimental demonstration would seem 
desirable. However, a practical difficulty is that real rotors ex- 
change energy with their surroundings for a variety of reasons: 
radiant heating or cooling, transfer of momentum to gas molecules 
exerting a frictional torque on the rotor, interaction of the rotor 
with external electric or magnetic fields, and so forth. 

Neglecting astronomical objects, an evaluation is made of the 
feasibility of reducing known torques to a level where unknown 
torques may be detected. One unknown torque considered here, is a 
gravitational radiation by the rotor, causing an energy loss or ro- 
tational decay by the rotor. With present instrumentation, certain 
kinds of gravitational radiations now seem on the verge of detecta- 
bility. 


INTRODUCTION 


Save for the frictionless planets, the concept of a free rotor was 
long considered academic. Real rotors always suffered the over- 
riding effects of bearing friction and chatter. The introduction of 
the air bearing was a vast improvement, although gas friction still 
drags heavily on the rotor. During the last twenty-five years, the 
work of J. W. Beams and his students (1) has produced a technique 
for suspending a small rotor magnetically in a vacuum. (See Fig- 
ure 1.) Under ideal conditions, this method provides a frictionless 
mounting, and small steel balls and cylinders can be spun magneti- 
cally to the bursting speed. 


This work was supported in part by U.S. Army Transportation 


Research Command, Fort Eustis, Virginia. 
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Upon turning off the drive motor (rotating B-field) the rotor will 
coast freely for many days and transfer its kinectic energy to its 
surroundings only gradually. It is the manner in which a small steel 
sphere or cylinder loses its energy that concerns us here. 

The kinectic energy of a rotor is of course 

Se : 

E=5 lw = 5 I(2nf) (1) 
where I is the moment of inertia and w and f are the frequencies 
expressed as radians per second and revolutions per second re- 
spectively. The power loss of the rotor is then 











dE/dt = Iw dw/dt = (27)*I df/dt (2) 
and 
4, y/ 
dE/dt _ 9 dw/dt _ 2 df/ dt (3) 
E w f 
The factor we — or os has the units of reciprocal seconds 


and is called the “decay ratio.” It is a measure of the rotor’s de- 
celeration at any frequency and is used throughout this paper, 
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because of its convenient interchangeability of w and f. Note that the 
decay ratio is independent of the shape, mass and hence, moment of 
inertia of the rotor. 


I Gravitational Radiation 

Table 1 outlines the various kinds of g-radiation. The first type 
(dg/dt # 0) or Cavendish (2) radiation is caused by a time variation 
of the g-force and is most often characterized as the quadrupole ra- 
diation emitted by a rotating dumb-bell. The second type (dg/dt = 0) 
is hard to visualize but always seems to be some function of the 
(v/c) ratio. 


Table 1 - Gravitational Radiations 


TOTAL POSSIBLE POWER 


_ 
ef | RADIATED 
FIRST TYPE (Quadrupole) kw 
SECOND TYPE 
LAPLACE (Doppler Shift) k’( ~ ) Ew 
BIRKHOFF (Accelerational) k"(=)° Ew 
EINSTEIN (Eddington) k"(—)° Ew 
GENERAL FORM k(~) Ew 


Laplace radiation (v/c) is a type of gravitational Doppler effect 
or aberration drag (3). Birkhoff radiation (v/c)* is conceived to 
arise because each point mass acts in such a way as to detract from 
the energy it possesses by virtue of being in an acceleration field 
(4)(5). Einstein radiation (v/c)* has been discussed by Eddington (6) 
who derived the magnitude of its effect on rotating cylinders. 

Using Birkhoff radiation, the (v/c)* type, as an example, we may 
evaluate the decay ratio, s , for a cylindrical (k” = 4.57) or 

Ww 
spherical rotor (k” = 3.68). Since Beams has shown (7) that steel 
rotors burst at a maximum rim velocity, v, of about 10° cm/sec, 
regardless of size, we may quickly estimate the maximum g- 
radiation effect for a cylindrical rotor. 
dw/dt _1,,,v,° 
= k (=) w 


oy 
(3 x 10! °y 
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df/dt 


«16 
7— = 5.3x10 (f) (5) 


The smallness of this power loss can be emphasized by imagining a 
megacycle rotor (10° revolutions per second) with a decay ratio of 
the order of 107° (sec™’). This means the rotor would lose 107* 
rev/sec, every second, or would require 10* seconds (about 3 hours) 
to lose one revolution per second. 

When the decay ratio tor g-radiation is compared to the decay 
ratio due to other causes (gas friction, etc.), the task of experimen- 
tally discriminating out the g-radiation effect seems formidable 
though not impossible. 

It might be well to consider the Laplace type since this would 
give a more favorable v/c factor than Birkhoff radiation. However, 
long standing astronomical observations, which have never been 
challenged, give little support to this aberration drag. For example, 
it was Laplace himself who showed that observable effects should 
occur in the motion of the moon-of a kind that are in fact not ob- 
served - unless c, the speed of gravitation, was at least 100 million 
times the speed of light (3). 


II. Uncertainty in Frequency Measurements 








A. Uncertainty Principle 


Small changes in rotational frequency might be significant or 
might be caused by random errors of measurement. In order to 
rely on frequency changes, it is necessary to explore the precision 
of frequency measurement and the limits of data refinement. 

The ultimate, inherent limitation in any measurement is, of 
course, the Heisenberg Uncertainty Principle, which leads directly 
to the concept of quantum states. An object must have a statistical 
number of quantum levels available to it, before its frequency can 
be measured with any degree of precision. In the case of a solid 
cylinder rotating at high speeds, this means the angular momentum 
should not be allowed to become quantized. 

This means that as the radius of the rotor decreases, the fre- 
quency must exceed an ever increasing number to obtain a statisti- 
cal number of quantum states. Estimations quickly show that the 
effect will not be serious for rotors larger than 10° cm in diameter. 


B. Effect of Light on the Rotor 


As a consequence of the Uncertainty Principle, it must be re- 
membered that any attempt to measure rotational velocity too 
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exactly would result in a disturbance of this velocity by the measur- 
ing probe and method of measurement. Presumably, the least dis- 
turbance would be caused by a light beam, when the rotational fre- 
quency is measured photo-electrically. This disturbance might not 
be negligible, and it might be well to estimate the effect of light on 
the rotor. 

The effect of light may be subdivided into absorption, trans- 
mission (absorption with re-emission), and reflection. These first 
two effects would disturb the temperature of the rotor, while the 
last would cause a pressure On the side of the rotor. Let us con- 
sider the effect of reflection, first. 


1. Light Pressure. In order to measure the frequency of ro- 
tation, the photomultiplier tube must receive at least one photon of 
reflected light from the rotor per revolution, If the efficiency of the 
tube is one per cent, this means that 100 photons per revolution 
must be received as a bare minimum. To be safe let us assume that 
the rotor reflects 1000 photons per revolution. As can be seen from 
Table 2, this could lead to a fairly intense beam for high frequency 
rotors. 


Table 2 - Light Pressure 





Frequency Light Reflected Force Pressure 

(rev/sec) (photons/sec) (ergs/sec) (dynes) (dynes/cm* )(mm Hg) 

40° 10° 4x10~* 8/3x107* 3x10" 2x10" 
10° 10° 4x107* 8/3x10" 3x107*° 2x10" 
10” 10*° 4x107* 8/3x10" 3x107™' 2x10~*° 
10° 104 4x10~" 8/3x10""  3x107™* 2x107" 
10° 10? 4 8/3x10".,  3x10™* ~~ 2x10™* 


It may be seen from Table 2 that this pressure is negligible com- 
pared to the vacuum pressure for ordinary rotational frequencies. 
However, above 10 revolutions/sec, light will begin to exert sig- 
nificant pressure and/or friction on the rotor. In other words, we 
will not be able to measure the rotational velocity without disturbing 
this velocity. The recovery time of an ordinary photomultiplier tube 
is about 107’ sec, and this limits the frequency to 10° rev/sec any- 
way. 

2. Radiant Heating. If the rotor is something less than a per- 
fect reflector, it will absorb some of the incident light and become 
warmer. As it becomes warmer, the rotor will begin to radiate 
more and more energy to the cooler surroundings. This radiation 
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follows the classical T* law. The emissivity coefficient, e, may 
vary from 0.05 for a very bright surface (such as an iron whisker) 
to a value of 0.5 for a dull, cast iron finish. 

It should be obvious that the rotor will heat up to the tempera- 
ture where the rate at which it emits energy will approach the rate 
at which it absorbs energy. This radiative equilibrium temperature 
may be calculated by using the assumed light reflection values in 
Table 2. The values for light absorbed are listed in Table 3 and 
calculated from the relation (6). It should be remembered that 
higher frequencies require smaller rotors to prevent bursting from 
centrifugal strain. 

Light absorbed = (7— 
Table 3 - Radiant Heating 
Frequency Light Reflected Light Absorbed Area T of Rotor® 


(rev/sec) _ (ergs/sec) _(e=0.5) (e=0.05) (cm* ) (e=0.5) (e=0.05) 


10° 4x10™" 4x107* 2x10°” 4x107' Small Small 
10° 4x10" 4x10°° 2x10°* 4x107* small small 
10" 4x10" 4x10"* 2x10°°  4x10™° 0.32°C 0.17°C 
10° 4x10" 4x10" 2x10" 4x10"" 156°C 101°C 
10° 4 4 2x10 4x10" large large 


) 
‘The temperature of the surroundings is taken at 27 C, 


As can be seen, the rotor will not be heated appreciably until 
rotational frequencies in excess of 10° rev/sec are sought. The 
rotor will react to this heating by changing its dimensions or 
moment of inertia, and hence, its rotational velocity. Once again the 
measuring instrument will disturb the thing it is measuring. 

The exact effect of temperature changes on rotational velocity 
will be discussed ina latter section. For the time being, it might 
be well to note that 10’ rps seems to bea practicable limit. For 
frequencies below this limit, the measuring light beam should be 
left on continuously so that the rotor may be brought to radiative 
equilibrium before the critical measurements begin. 


ITI, Gas Friction 


It is anticipated that the magnetically suspended rotor will be 
spun in a vacuum of 107° to 107*° mm Hg. However, even at these 
very low pressures there may be sufficient gas to cause a signifi- 
cant frictional drag on the rotor. 


) (light reflected) (6) | 
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When dealing with a high vacuum, it must be remembered that 
the mean free path of the gas molecules is large. The normal vis- 
cosity coefficient, characterizing internal friction of the gas, does 
not apply. Instead a free-molecule viscosity coefficient, Z, is em- 
ployed (8). The value of Z is given by 


M.W. 3 


Z = P(oTRT”” (7) 


where P is the gas pressure; M.W. is the molecular weight; T is the 
absolute temperature; R is the gas constant. Note that Z, P, M.W., 
and R are all in c.g.s. units. 

The decay ratio for gas friction on a spherical rotor is then 

dw/dt _52Z (8) 
w rp 
where P is the density of the rotor. Equation (8) has been verified 
experimentally by Beams (7). 

A range of values for the decay ratio from equation (8) are given 
in Table 4. It is assumed that a steel sphere ( p = 7.8) is rotating in 
air at very low pressures and 300° Absolute. By comparing the de- 
cay ratio caused by gas friction to that caused by other factors, it 
will be seen later that gas pressures will have to approach 107° 
mm Hg before frictional drag becomes negligible. 


Table 4 - Decay Ratio Due to Gas Friction on Spherical Rotors 





Pressure |__Radius 1 107" = 10° : ; 

(mm Hg) | r(cm) aS OD, =. 4) Baa e - 
107° 1.16x10™" 1,16x10"° 1,16x10" 1.16x10"~ 
107° -8 -7 -6 =5 

10 10 10 10 

107" 107 - 107° 107. 107° 
107° 107 107° 107° 107" 
107° 10". 107*° 107° 107° 
107° 10” 107" ” is 107° 





IV. Rotation of the Earth 
The following axioms may be observed: 
1. The axis of spin of the rotor (gyro) tends to be vertical due 
to gravitational attraction, It seems the point of support of sucha 
rotor would logically be above the center of gravity, since the lines 
of magnetic flux diverge from the bottom of the supporting solenoid. 
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2. Assume the direction of spin is counterclockwise looking 
down from above the rotor. 

3. An applied external moment cannot change the rate of spin of 
the rotor or flywheel; it increases the angular momentum of the 
mounted wheel and ring by making it precess about the vertical] 
axis (9). 

4. Aspinning gyro tends to maintain its direction fixed in space, 
Hence, the rotation of the earth produces a forced precession about 
the precession axis, which is a vector parallel to the axis of rotation 
of the earth passing through the center of support of the rotor. 

5. This enforced precession calls forth a gyroscopic reaction 
moment that is equal and opposite to the external moment that would 
be required to produce the forced precession. 


The rotation of the earth influences rotational decay only in an 
indirect fashion. A forced precession is elicited in the rotor, which 
adds or subtracts less than one revolution per day on the absolute 
rotational velocity. However, once the rotor has attained 2 maxi- 
mum spin, it is only the time rate of change of this spin that will be 
measured, and the earth’s rotation will be cancelled out. 

The force exerted on the rotor by the support coil, causing it to 
precess, will not be reflected by a gain or loss in spin. Rather this 
work done by the coil will be transformed into the kinectic energy of 
precession, while the spin remains constant. 

The earth’s rotation does, however, induce a magnetic drag on 
the rotor. This is because the axis of spin is inclined to the axis of 
magnetic support due to a gyroscopic reaction moment. The effect 
of this and other horizontal components of magnetic fields will be 
considered in the following section. 


V. Estimation of Magnetic Drag 





By total energy considerations, one may estimate the drag on a 
cylindrical rotor (radius, r, and length, 1) caused by a horizontal 
magnetic field. 


€=vXB =rwcos¢?B (9) 
(mks) x 

where ¢ is the electric field strength (volts/meter); v is the rim 

velocity (m/sec); B is the magnetic field strength and B, is its 

horizontal component (webers/m’*); where @ is the angle of revo- 

lution; where dé/dt is w. If J is taken as the current density 

(amps/m’*) and o as the conductivity (1/ohm-m), then 
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= 

Power loss = | dv = { oe?dv 

A CG 
= { f[ o(rw)? cos? ¢ B? rd@ dr 7 (10) 

» 9 

dE/dt = 2 dw, dt » = B*. 
E Ww p 


Using appropriate values for the density (kg/m*) and conductivity of 
steel and converting webers to gauss 


dw/ dt 


A 


= 4.27x10"°B?2 — (gauss) (11) 


The significance of equation (10) may be realized if it is re- 
membered that the horizontal component of the earth’s magnetic 
field is about 0.2 gauss. Then too, the earth’s rotation, as we have 
seen in the last section, may indirectly cause a horizontal magnetic 
field component of as much as 3 gauss across the rotor. According 
to these estimates then, magnetic drag is far from negligible and 
will have to be carefully considered. 


VI. Changes in Temperature 


Since metallic rotors will expand or contract with changes in 
temperature, the law of the conservation of angular momentum 
would require that this expansion or contraction be accompanied by 
corresponding changes in angular velocity. If rotational frequencies 
can be measured to one part in 10° or 10’, then a change in tem- 
perature of 1° C will make a measurable change in frequency. If a 
rotor is heated up to 100° above room temperature while being 
driven up to 10° rev/sec by the magnetic drive coils, it will gain 
240 rev/sec upon cooling to room temperature. Then too, it can be 
seen from Table 3 that radiant heating by a light beam may slow a 
small rotor down considerably. 

These considerations are not of overwhelming importance be- 
cause of the lack of a heat transfer variable. For example, upon 
cooling from a higher temperature, the acceleration of the rotor can 
not be calculated because the cooling rate has not been specified. 
What is needed is some type of acceleration factor, dw/dt. In order 
to evaluate this factor, we must know the heating or cooling rate. If 
heating or cooling is achieved by radiation only, T* law may be used 
to calculate the energy transfer rate. Table 5 gives a range of 
decay ratios due to thermal fluctuations for both cylindrical and 
spherical rotors of various diameters. 
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Table 5 - Effect of a Thermal Fluctuation 
of + 1° C on Decay Ratio 


Diameter Spherical Rotors Cylindrical* 


(cm) (e=0.5) (e=0.05) (e=0.5) (e=0.05) 

2 6.2x107* 6.2x107*° 6.5x107° 6.5x107~” 
2x107° 6.2x107° 6.2x107° 6.5x107*° 6.5x107° 
2x107? 6.2x107" 6.2x107° 6.5x107" 6.5x107° 
2x107° 6.2x107° 6.2x10~" 6.5x107° 6.5x107~" 
2x107* 6.2x10~° 6.2x107° 6.5x107° 6.5x107° 
A mmstant Snape was issumeada tor the <« iinaer where 
(7 = 4 


VIL Elastic Strain 


For an estimate of this hoop stress, we may follow Beams and 
take it approximately equal to r°w* or N*. The change in radius 
may then be calculated from the change in circumference of the 
rotor, using Q as the circumferential tension. 


Young’s Modulus (steel) = 2x10** (dynes/cm’ ) 
= Q/(A(27r)/2rr |} (12) 


The significance of a variable radius for the rotor is this. Im- 
plicitly contained in all that has gone before was the assumption that 
any power loss of the rotor was accompanied by a decrease in ro- 
tational velocity while the radius remained constant. That is 

dE/dt/E = 2 dw/dt/w (13) 
where E for a cylinder is  mr* w’. 
Now, however, we must write 
dE/dt/E = 2 dw/dt/w + 2 dr/dt/r (14) 
since r is also a variable. The problem is to estimate what contri- 


bution the r-term in equation (14) will make to the total power loss. 
From equation (12) we may write 


Ar/r = Beto = r? wy” /2x10' (15) 


where r, is the radius of the rotor at rest. Then 
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. w? 
Su P= st 107 [3 dr/dt/r + 2 dw/dt/w 


a6 Cx 
2 rw (16) 
dr/dt/r _ 2 x 10 ea. oe 
dw/dt/w 1. 3r 10** - 1.5 r*w* 


2x 10” 
Equation (16), then, gives us the relative contribution of the change 
in radius of the rotor compared to the change in angular velocity. 


A range of values are given in Table 6. 


Table 6 - Relative Importance of Change in Radius 


Radius 1 107? 107? 107° 107* 

_ (cm) 

Frequency 

(rev/sec) 
105 9.7x107' 3.9x10™* 3,.9x10"° 3.9x107"  3.9x107° 
10° 9.7x10°* 3.9x10™* 3.9x10"*  3.9x107” 
10” 9.7x107' 3.9x10"* 3.9x107° 
10° 9.7x107* 3.9x107* 
10° 9.7x107* 


The values recorded are the right hand side of equation 16 and 
hence, represent the fractional contribution of the change in radius 
to the power loss. It can be seen that in most cases onlya negligible 
amount of the rotor’s energy is stored in mechanical strain. In 
these cases, we are justified in taking the decay ratio (dw/dt/w) as 
a measure of the fractional power loss (dE/dt/E). However, for 
larger rotors at higher frequencies, mechanical strain plays an in- 
creasingly important part until it eventually comes to equal terms 
with angular velocity. The blank spaces in Table 6 indicate rim 
velocities which are well above the bursting strength of known ma- 
teriais. 

It would seem then that a rotational frequency of 10° (rev/sec) 
once again represents a natural limit, since for higher values the 
rotors become vanishingly small. Not only that, if mechanical strain 
starts to play an important part in power loss, the rotors velocity, 
shape, and dimensions will have to be measured simultaneously. It 
would probably be well to do this anyway and be prepared to explain 
anomalous results. A bulging of rotor at the center would give dif- 
ferent results than those predicted here. 
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VIII. Conclusions 


From what has been said above, it may be seen that some of the 
factors considered may be vital while others may be ignored. Of 
primary concern isthe temperature control of the rotor’s surround- 
ings and the production of an ultra-high vacuum in which to spin the 
rotor. The temperature of the chamber should be controlled +0.01'C, 
and gas pressure should be below 1x10~*° mm Hg. Both of these 
conditions are on the fringe of what is now possible to produce ina 
laboratory. 

Factors which are no less critical but more easily nullified are 
magnetic drag, rotation of the earth, and mechanical vibrations. The 
last two factors result ina form of magnetic drag. However, there 
seems to be a good chance of dampening vibrations and controlling 
stray magnetic fields with a Helmholz coil, Then too, it should be 
noted that the decay ratio due to g-radiation is a function of angular 
velocity whereas the magnetic drag decay ratio is not. By going to 
smaller rotors with higher permissible frequencies, the magnetic 
drag effect is made less dominant. 

Factors which may be controlled by choosing the proper shape, 
size, strength, and surface condition of the rotor are radiant heating 
and the stress energy of the rotor. Both of these factors seem to 
place an upper limit on the rotational frequency of one hundred 
megacycles. Table 7 gives a comparison of the decay ratios ex- 
pected from the various factors discussed here. 


Table 7 


Factors Influencing Megacycle Rotor 


r 7 . T T . 


Factor Condition Decay Ratio 
Se : = . = 
[o- Radiation 10° rps atl a 
Gas Friction 107*°mm He isi" 
Temperature of +0.1°C. 3x 107° 
Surroundings 
| Light Pressure 1000 photons/ rev ~107* mm Hg 
RFE: | Barth’s Field or Tilt ee: 
| Magnetic Drag ef Mater 10 


=m ——+t - a= p= 
| Mechanical Strain | Steel at 10° rps <1% error 


Merckens et one aA ocaetact — 


Notice that the effect of temperature is cyclic and will cancel 
out over long periods of time. Hence, even a temperature fluctuation 
of 71.0°C could be tolerated if the rotor is observed over a period 
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of 24 hours or so, [If the light beam, heats the rotor appreciably, it 
can be left on continuously, as mentioned above, so that the rotor 
can come to radiative equilibrium with its surroundings. 

However, the one irreducible problem that remains is the tilt of 
the rotor. As explained above, the distance between the magnetic 
center of support and the center of gravity of the rotor is not known, 
In lieu of performing the experiment at the north or south pole, a 
laboratory investigation of center of support and magnetic drag has 
been inaugurated here and is presently in progress. 
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pp 37-49, Interscience Publishers, New York, N. Y. (1958). 
(10) The exact angle of tilt due to the gyroscopic reaction moment 
of the rotor is difficult to ascertain because the lever arm 
must be known, which is the distance between the center of 
support and the center of gravity of the rotor. Assuming an 
extreme case of the center of support being at the top of the 

rotor, we get 

Sin — = (3.19 X 107°) rw = 3.19 x 107° 
where rw is the rim velocity (10°cm/sec) and — is the small 
angle of tilt of the rotors axis from the vertical. Multiplying 
the support field strength (which could be of the order of 
10° gauss) by Sin & gives the horizontal field component. 
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A “Chain Block” Experimental Design with 
An Application in Fuel Testing * 


By Kan-chen Peng 
Research Laboratories 
General Motors Corporation 


1. INTRODUCTION 


An experimenter in the physical sciences often encounters 
situations in which 


(1) The number of “treatments” (or tests) is large and may 
considerably exceed the size of “block” (such as the length 
of a day, the capacity of a certain set of test equipment, etc.). 

(2) Within the blocks, comparisons are relatively precise. 

(3) Elimination of effect due to blocks is desired. 


The classical experimental design may fail to meet these con- 
ditions successfully. With these needs in mind, Youden and Connor 
[1] developed the simple chain block designs. Mandel [2] has shown 
that some of these designs can be generalized so as to adjust for 
two-way (columns and rows) heterogeneity. This paper discusses 
the essentials of these designs with an application in fuel testing, 
and indicates that these designs can be used together with a multiple 
comparison procedure [3]. 


2. DESCRIPTION OF THE DESIGN 


The principle of the simple chain block design may be illustrated 
by an example in Table I, for 17 treatments in blocks of 6 or 5 units. 

The letters denote the treatments. Capital letters refer to treat- 
ments having two replications; small letters refer to treatments 
with a single replication. 

The design has a basic linking property. C and D appear both 
in Blocks 1 and 2, and form the link between these two blocks. E and 
F, Gand H, I and J, and A and B form the links between Blocks 2 


*The author is indebted to many people at the Research Labora- 
tories of General Motors Corporation, especially to Dr, F. W. 
Bowditch for proposing the problem that led to this paper, and to 


Mr, D. E. Hart and Mr. J. E., Dallemand for encouragement and 


helpful discussions. 








and 3, 3 and 4, 4 and 5, and 5 and 1 respectively. Therefore, the 
design may be regarded as a chain. The linking property is used to 
compute the block effects. Each treatment mean is then adjusted 
for the effects of the block or blocks in which it lies. 


TABLE I 


A Simple Chain Design for 17 Treatments 
in Blocks of 6 or 5 





BLOCK 
l 2 3 4 5 
A Cc E eit 
B D F H | J 
Cc’ E’ G’ I A’ | 
D’ F’ H’ J’ B’ | 
k l m n 9) | 
P q | 




















The design is flexible. The basic groups such as A, B may have 
1, 2, 3, or more letters (or treatments). However, all groups should 
have the same number. Additional treatments with a single replica- 
tion like k, 1, m, may be included, subject only to the limitations 
imposed by the size of block, 

The number of error degrees of freedom for a design is 


N-v-b+l 


where N is the total number of letters counting replicates, yp the 
number of letters ignoring replicates, and b the number of blocks. 
The design in Table I has 6 error degrees of freedom. This number 
should be checked before adopting a design. 

Mandel [2] has shown that some of the simple chain block de- 
signs can be generalized in such a fashion that elimination of bias is 
achieved not only for blocks (columns), but also for another factor, 
which may, for example, be identified with order within blocks 
(rows). Table II gives a design with this property for 8 treatments. 
If the columns are regarded as blocks, this is a typical chain block 
design, and the chain goes 1, 3, 2, 4, 1. 

Again, if the rows are regarded as blocks, the chain goes I, IV, 
II, Ill, I. Thus the arrangement is a chain block design for both 
rows and columns, and treatment means may be corrected to re- 
move the effects of the columns and rows. 











Now 
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TABLE II 
A Generalized Chain Block Design 


ROWS 1 2 
. » | 

E | 

o lied 

c' | p’ | 

L 

Other similar designs are given in Mandel’s paper. Note that 


the number of blocks must be even and that the number of treat- 
ments must be a multiple of the number of blocks. 





mwa ale 
mr xt Ole 


IV 


——-----———_ 











3. BASIC METHOD OF ANALYSIS 


An outline of the method of analysis, needed in Section 4, is 
given here. This follows the presentation developed in [1] and [2]. 

Let a. ,a pM, neers a. denote v treatments, or more generally, 
groups of treatments with each group having the same number of 
treatments. Let a; denote the average yield of group aj of the first 
replication, and a; denote the average yield of group a; of the second 
replication. Then, the general form of a simple chain block design 
may be shown in Table III. Note that this corresponds to the upper 
part of Table I. 





TABLE III 
The General Form of a Simple Chain Block 
— a a ne eS 
““._-—s«;BLOCK 

1 2 3 i V 
REPLICATION ~_| 

I a a, ee ere 

II a! By Ba wees Bigs os0 B, 

Now, define the following relations 
(1) P; =a',, - ay, (j+1 is interpreted mod v) 
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Vv Vv 
(3) G = grand sum = ai + Sai 


i=1 i=1 


(4) Sequence S = (v-1), (v-3), (v-5), ...., -(v-5), -(v-3), -(v-1) 
where S represents an arithmetic progression with common 
difference -2, 


Let e; (j = 1, 2,...., v) represent the bias of block j, that is, e; is 
the systematic error affecting each yield in block j. If all system- 
atic errors are taken with reference to the overall average, we may 
impose the restriction, 
Vv 
(5) >> e; =0 


Then, to compute 6, which is the estimate of an ex, we permute the 
Py. cyclically such that P, is the first term, and write down the P’s 
below S: 


S sequence v-1_ v-3 v-5 eee -(v-5) 0 -(v-3) )~—s -(v-1) 
Pvalues P, Pri Pris nooo Mice Py.-2 Pinas 


Now, form the sum of products B, as follows: 

(6) By = (v-1)(Py, - Py-,) + (v-3)(Pu, - Px-2) 

+ (v-5)(Prio- Pry) +... 

where the subscripts are interpreted mod v, and finally, 

(7) 6, = : B 

k = oy ~ 

The recursion formula 

(8) @x4,=@y- Py +P 


which can be obtained from (6) and (7) is useful for computing é,, 
é,, ..., after 6, is computed from (6) and (7). 

The analysis of variance is shown in Table IV, as given in [1] 
and [2], where derivations of least squares solution of the analysis 
can be found. The analysis for the generalized chain block design 
follows immediately from the discussions in Section 2 and the for- 
mulas given above. Only proper rearrangement of the columns and 
rows is required. The computed biases for columns and rows, with 
signs changed, are applied to the original observations as additive 
corrections. 
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TABLE IV 


Analysis of Variance of Simple Chain Block 








Degrees 
Source of Sum of Squares 
Freedom 





n- Total 2v -1 S, = 
Treatment Ignoring Blocks v-1 S, = 


Blocks Eliminating Treatments v-1l S, = 











he Rows &, = = 
sg ov 








4. AN APPLICATION IN FUEL TESTING 


A generalized chain block design and a multiple comparison 
procedure have been successfully applied to a fuel testing problem 
which was proposed by Fred Bowditch, Fuels and Lubricants De- 
partment, General Motors Research. 

Briefly, the problem involves testing eight different fuels ona 
single test engine. Experimental conditions made it necessary to 
extend the tests over discrete test periods during a day and several 
days, and it was desired to remove the bias from days and periods. 
A design as shown in Table II was used for this purpose. A set of 
16 measurements from a total of 64 measurements is shown in 
Table V. These measurements represent crankangle degrees before 
top dead center (BTC). 


TABLE V 


Fuel Testing Measurements (° BTC) 
in a Generalized Chain Block 


DAYS 





_— 
I (A) 100.78 (B) 97.98 (C) 99.10 (D) 101.78 
II (E) 105.52 (F) 88.55 (G) 90.95 (H) 109.00 


Ill (G’) 75.08 (H’) 110.65 (B’) 112.58 (A’) 114.92 





moor a Dv 


(D’) 103.60 (F’) 88.25 (E’) 115.92 
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After the effects of columns and rows are removed, the differ- 
ences between the two replications of all fuels are considerably 
reduced as shown in Table VI. 
































TABLE VI 
Peet Uncorrected | |4| Corrected | |4 
~A | 100,78 | 114.92 | 14.14 | 113.67 | 110.88 2.79 
B | 97.98 | 112.58 | 14.60 | 100.57 | 103.36 | 2.79 
Cc | 99.10 | 87.48 | 11.62 | 96.39 97.41 | 1.02 
D | 101.78 | 103.60 | 1.82 | 104.25 | 103,23 | 1,02 
E | 105.52 | 115.92 | 10.40 | 112.64 | 115.43 | 2.79 
F 88.55 | 88,25 30 | 85.37 82.58 | 2.79 
G 90.95 | 75.08 | 15.87 | 82.47 | 81.46 | 1.01 
H 109.00 | 110.65 | 1.65 | 105.70 | 106.73 | 1.03 
TABLE VI 
Fuel Testing Measurements 
DAYS 

‘ . 4 Me) sl) 

1} (A) (B) (c) (vo) | (A,B) | CD) 

| 100.78 | 97.98 | 99.10 | 101.78 | 99.38 | 100.44 

a}. (F) | (G) (H) | (E,F) | (GH) 

105.52 | 88.55 | 90.95 | 109.00 | 97.04 | 99.98 
= 7 oe | 1‘ ia 119.98 11499 “oaee T1878 | 
PERIODS ; ' | 
IV (C’) (D’) (F’) (E’) |(C’,D’)|(E Ay 
87.48 | 103.60 | 88.25 | 115.92 | 95.54 | 102.08 | 








I+1I | (A,E) | (B,F) | (C,G) | (D,H) 
2 | 103.15 | 93.26 | 95.02 | 105.39 


IlI+IV |(C’,G’)|(D’,H')|(B’,F’)|(A’,E’) 
2 | 81.28 | 107.12 | 100.42 | 115.42 
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This suggests that the design in this case was effective in re- 
moving biases. This conclusion is confirmed by the analysis of 
variance. The detailed procedures for removing biases are shown 
in Tables VII to X, and the analysis of variance in Table XI. The 
computations are based on the notations and formulas in Section 3, 
with rearrangement of rows and columns as indicated in Table VIII 
and Table IX. 

A multiple comparison procedure based on Duncan’s method [3] 
is then used on the 64 corrected measurements to rank the fuels. 
Note that 6 degrees of freedom are used to correct biases for each 
of 4 sets of 16 measurements, This fact has to be considered before 
the multiple comparisons can be made. The application of Duncan’s 
method is rather straightforward. The reader is referred to Dun- 
can’s Original paper for detailed information of the method. 


TABLE VIII 


Column Effects 
































1 3 2 4 
= 3 2 3’ 4’ 
[+--+ 
| (A,E) (C,G) (B, F) (D,H) 
| 103.15 95.02 93.26 105.39 
| €’.c’) (B’,F’) (D’ .H’) (A’,E’) | 
| 81.28 100,42 107.12 115.42 
Pell Miaseneatt Tinsinandt ——_ 
P \(C’,G’-C,G) | (B’.F -B, F)|(D’ ,H’ -D,H) (A’ ,.E’-A,E) | 
| -13.74 7.16 1.73 12.27 |P = 1.8550 
5 3 1 -1 -3 | 
<a — alieeeeidaall 


= -9,0750 - (-13.74) + 1.8550 = 6.5200 
Yy = ¥,= 6.5200 - 7.16 + 1.8550 = 1.2150 
y 


~ 
nN 
' 
~> 
wo 
\ 


= 1.2150 - 1.73 + 1.8550 = 1.3400 


Ye 
(1, is computed by (7) in Section 3, and the other ¥i,’s are com- 
puted by (8).) 








? 


? 
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TABLE [Ix 


Row Effects 
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I Ill II IV 
j 1’ 2’ 3’ 4’ 
(C,D) (A’ ,B’) (G,H) (E',F’) 
100.44 113.75 99.98 102.08 
(A,B) (G’ ,H’) (E, F) (C’,D’) 
99.38 92.86 97.04 95.54 
P |(A,B-A’,B’) | (G’,H’-G,H)|(E,F-E’,F’)|(C’,D’ -C,D) & 
-14.37 -7.12 -5.04 -4,90 P = -7,8575 
Ss 3 1 -1 -3 
p,, =p, = ; [3(-14.37) + 1(-7.12) - 1(-5.04) - 3(-4.90)] = -3.8112 
P., = Pyy, = -3.8112 - (-14.37) - 7.8575 = 2.7013 
P., = Py; = 2.7013 - (-7.12) - 7.8575 = 1.9638 
P,» = Pry = 1.9638 - (-5.04) - 7.8575 = -.8537 
p,’ is computed by (7) in Section 3, and the other Y; ’s are com- 
puted by (8).) 
TABLE X 
Calculation of Corrected Values 
1 2 3 4 1 2 3 4 
Correction | +9.08} -1.22) -6.52] -1.34 
I +3.81/100.78| 97.98) 99.10/101.781/113.67/100.57| 96.39 | 104.25 
II ~1.96}105.52| 88.55] 90.95 /|109.00}|112.64| 85.37] 82.47 |105.70 
Il -2.70| 75.08) 110.65 |112.58 |114.92]| 81.46/ 106.73 | 103.36 |110.88 
IV + .85| 87.48! 103.60| 88.25 |115.92|| 97.41|103.23| 82.58 |115.43 | 
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TABLE XI 
Analysis of Variance of Fuel Tests 
- Degrees of | Sum of Mean no 
= panes Freedom Squares Square y-matie 
ae 
a Total 15 2009. 3676 -- -- 
Fuels 7 1552.234 221.7477 25.117 
Columns 3 379.8350 126.6117 14,341 
Rows 3 59.6417 19.8806 2.252 
— Error 2 17.6569 | 8.8284 -- 
i75,) The sum of squares for Columns is obtained as follows. 
i 2 {; |(-9.0750)(-13.74 - 12,27) + 6.5200(7.16 + 13.74) 


— + 1.2150(1.73 - 7.16) + 1.3400(12.27 - 1.73)]} 
= [236.0408 + 136.2680 - 6.5974 + 14,1236] = 379.8350 
The sum of squares for Rows is obtained as follows. 
2 {5 [(-3.8112)(-14.37 + 4.90) + 2.7013(-7.12 + 14.37) 
F + 1.9638(-5.04 + 7.12) - .8537(-4.90 + 5.04)] } 
= 36.0921 + 19.5844 + 4.0847 - .1195 = 59.6417 


(The factor 2 is used, here and above, because the data in Tables 
VIII and [IX are averages of two original observations.) 


The sum of squares for Total is 


(1602, 14)° 
2x8 


162437.6536 - 160428.286 = 2009. 3676 


162437.6536 - 


4 The sum of squares for Fuels (uncorrected) is 


161980.520 - 160428.286 = 1552.234 
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A Measure of the Optical Quality of 
Proposed Automobile Windshield Designs 


By J. R. Gagnon 

Systems Engineering Department 

Engineering Research and Advanced Product Study Office 
Ford Motor Company 


ABSTRACT 


A measure of the optical quality of automobile windshield de- 
signs is obtained by studying the refracting characteristics of the 
proposed glass surfaces. The deviation of light rays converging to 
an observation point is calculated for various directions of obser- 
vation in order to determine deviation patterns that can be analyzed 
along significant scanning paths. 


INTRODUC TION 


Conventional optical techniques can be used to determine the 
refracting characteristics of glass surfaces. However, in studying 
the optical performance of automobile windshields, it is advanta- 
geous to eliminate the need for actual glass models, and work di- 
rectly with the proposed dimensions. To present the problem and 
the selected approach for its solution, a schematic of the optical 
system will be explained. 

Figure 1 shows a windshield in the rectangular system of axes 
X’ Y’ Z’ whose origin O’ is at a point of observation determined by 
the location of a passenger in the vehicle. From that point, the 
entire windshield area can be scanned, and the quality of the image 
seen depends on the particular glass section through which light 
rays from an object are transmitted. 

Consider an arbitrary plane (P), in the environment, containing 
a vertical linear object ABC. A light ray from point B will meet the 
outer surface of the windshield at P,’, be refracted to P, on the in- 
side surface, and refracted again to O’. As will be shown later, 
the direction of the ray changes as it passes from one medium to 
another (air to glass or vice versa), and the resultant effect is that 
the ray P, O’ appears to come from a point B’ in the plane (P). 





> 
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Vector BB’ is defined as the deviation, in plane (P), produced by the 
windshield, when the direction of observation is O’P . 

A similar treatment can be done for a sufficient number of rays 
coming from the object ABC, and the optical theory will show that 
the magnitude and direction of the deviation of each ray can be dif- 
ferent. Rays from points A and C could appear to come from A’ 
and C’ respectively, and the linear object ABC would be perceived 
by the observer at O’ as a curve A’B/C’. 

The departure of the image A’B’C’ from the actual shape and 
location of the object ABC is a measure of the optical quality of the 
windshield when observed along directions O’P;. The mathematical 
development which follows will lead to formulas expressing the de- 
viation vectors as a function of design parameters. 


ANALYTICAL DEVELOPMENT 


The determination of deviation vectors in a convenient arbitrary 
plane will be done by computing the direction of light rays in the 
various optical media present. The law of refraction states that 
when a ray of light is incident on a surface separating two contiguous 
homogeneous media, it is refracted in the plane of the incident ray 
and the normal to the surface at the point of intersection and its 
direction in that plane is given by 

A A 


HW, sin aN = yp, sin BN (1) 


where wy, is the index of refraction of the incident medium, 
HM. is the index of refraction of the refracting medium, 


is the angle between a unit vector a in the direction of 
the incident ray and a unit vector N in the direction of 
the normal, and, 


DN is the angle between a unit vector b in the direction of 
the refracted ray and the unit vector N. 


To trace a light ray as it passes from one medium to another, it 
is therefore necessary to know the direction of the incident ray and 
the direction of the normal to the surface of separation at the point 
of intersection. The former can be selected, and the latter com- 
puted if the equation of the surface is known. To obtain this equa- 
tion, the coordinates of surface points defining a small section are 
measured from the proposed windshield design, and used in a sur- 
face fitting program based on the method of least squares. The 
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GEOMETRY OF THE RAY-TRACING PROCESS 
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Figure 2 


degree of the equation, for any particular small section, depends on 
criteria related to the optical problem. For simplicity here, assume 
that the coordinates of a number of surface points yield an equation 
of the form 


2 2 
O(x,y,Z) = ao +ayX+Aa,y +ayX +asy +agxy-2= 0. (2) 


All of these points can now be taken, one at a time, and the following 
procedure repeated for each one. 
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Referring to Figure 2, a point N, on the outside surface §S, is 
selected. Let the measured coordinates of N, in the designer’s 
coordinate system be xy, yn. and zy. These coordinates satisfy 
the equation of the surface. 
Therefore, 


@ = aot A&Xy + A,yy + AsXy + ay¥N + 45Xyn Yn - 2y. = O. (3) 


The direction of the principal normal to S, at point N, is given by 


0P/AxN 0?/ayn 
COS Qn. = —>~—~-,, COS bn. = ————— 
IN} C N, ? iN} C N 
- (4) 
oY? OZN 


and cosy, = 


where cos ayn, COs By and cosy, are the direction cosines of 
the normal, and Cy_ is the square root of the sum of the squares of 


the partial derivatives of @ with respect to xy , yn, and w,, or 


Cn, = [ (= + 

OXN, dyYN, dZN, 
The principal normal to S, at N, meets the inside surface S, at a 
point P,. The segment P,N, is the given thickness t of the wind- 
shield and P,N, is the principal normal to S, at point P, because 
the surfaces are designed parallel. 


Let the coordinates of P, be xp , yp and Zp. 
1 1 1 





XN - Xp 
By definition, cosa, = : , } 
YN, ~ YP, 
cos By ~ ; . > (6) 
ZN, - ZP, 
and, cos N, = ; J 
So Xp, = Xn, - t cosan , 
Yp, = yn, - t cos By, (7) 


and, Zp = Zn - t cosy. 
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Select a point of observation O’ of given coordinates x,, y, and Z, in 
the designer’s coordinate system, and let O’ be the origin of a sec- 
ondary coordinate system whose axes are parallel to the respective 
designer’s axes. OP, is the direction of observation. 


—_ 1 
; 2 2 = 
O’P, = d, = [(xp - x,)° + (yp -y,)° + (zp, -2,)°]2 , (8) 
and, by definition, 
Xp. - X, 
cos a@, = ’ 
l d, 
Yp ~ os 
cos 8, = " , (9) 
i d, 
zp -2Z 
and, cos y, = 
1 ’ 
d, 


where cos a, , cos 8, and cosy, are the direction cosines of O’P.. 
It is also convenient to express the direction of O’P by a horizontal 
angle of scan A and an angle of elevation E. 


ZE =<(90°- y,), (10) 
or 
cos E = siny . (11) 
Also cos A cos E = cos @, , (12) 
cos a, 
therefore, cog. A 2 auqe 5 , (13) 
sin 7, 


The linear segment O’P, defines the direction of observation, 
but it can also be considered as a light ray which has been trans- 
mitted through the windshield after experiencing a refraction at the 
outside surface and one at the inside surface. To find the direction 
of the incident ray, it is necessary to determine the path followed by 
light in the glass and apply the law of refraction at both surfaces. It 
will be assumed here that the air and the glass are homogeneous 
media. 

The ray of light, in the glass, that intersects surface S, at point 
P, , and is then refracted in air to point O’, is in the plane of P, O’ 
and the normal N.N,. The angle between this ray and the normal 
is given by equation (1), Let this incident ray in glass be P/P, , 
where point PY is the intersection of the ray with surface S,. 
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By equation (1), 
> M a i r 


sin P’P N = —sgnOP nN, (14) 
11) glass i 
angle O’P N, can be expressed as a function of direction cosines, 
thus 
sib 


cos O'R N, = cosa, cosay, + cos8, cosBy + cosy, cosy, . (15) 


Equation(15) gives angle O’P, N, as a function of known parameters, 
and angle P,’ P,N, can then be calculated from equation (14). 

Point P,’ is in the plane determined by points N., P and O’, and 
the coordinates of P/’ will satisfy the equation of the plane. Let this 
equation be 


A .x+By+C,z2+1=0. (16) 


Substituting the coordinates of each known point, equation (16) gives 


A, Xn, + Byyn, + C, 2n, = -1, 
A, Xp, + B yp, + C, Zp, = -1, (17) 
and, A,x +B,y, +C,% =-1. 


Having the three equations (17), the unknowns A,, B, and C, can be 
calculated. 


It is now necessary to find the coordinates of point P’ - Xp, yp 
and Zp. 
ee a 
P, PY - d, = (xp = Xp, a +r (Yp = Yp, y° r (Zp- Zp y }2 , (18) 
Xp - Xp 
and, COS @p = — d . 
JP ~ YP. 
cos Bp = i, (19) 
d, 
Zp - Zp 
COS Yp = q A, 


where cos a, cos fp and cos yp are the direction cosines of B B’. 
The angle between B’P, and N,N, has been expressed by equation 
(14), and 
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ro, 
cos 4 PN, = COS @p cos an, 


(20) 
+ COS Bp cOS By + COS Yp COSyy . 
N, N, 
Equation (20) can be written 
d,-cos P,’P, N, = (xp - Xp ) cos ay, (21) 


a (Yp - yp, )cos By, ° (Z, at Zp. ) cos YN, ’ 
which is a first function of Xp» Vp and Zp: 


The coordinates of P,’ satisfy equation (16) because P,’ is in that 
plane, so 


A.Xp + Biyp+C,2p+1=0. (22) 


1 


Equation (22) is a second function of Xp, Yp and zp. 


Point P,’ is also on the outside surface of the windshield, and its 
coordinates will satisfy the equation of the surface, so 


@=a,+a,Xp+a,yp+a,xp+a, yp +4, Xp Yp - Zp = 0 (23) 


becomes the third function of Xp, Yp andZp. 


Having three equations (21), (22) and (23), the value of the coordi- 
nates X»,Yp and Zp can be determined, 


Note: When the surface equation is of the third, fourth or higher 
degree, a first approximation for xp, Yp and Zp is obtained by re- 
placing equation (23) by the equation of the tangent plane to the sur- 
face at point N,. Numerical methods are then used to obtain the 
exact solution for the coordinates. 


The direction of the incident light ray is computed by applying 
the law of refraction at point P,’ on surface S,. The direction co- 
sines of the principal normal to S, at point B’ are given by 


09/8x p ) 
COS Gy, a —, | 
0¢/3yp (24) 
08 hy §*——— 
cos aN. S , 
09/azp 


and, COS yy =——-—, J 
2 
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where 








y*}2 , (25) 


angle N/P,’R between the ray PP, and the normal at P’ is given by 
— 
cos N;/P,'R = cos @ cos ay, 
(26) 


+ COS bp cos BN, + COS yp COS YN, » 


and the angle between the same normal and the incident ray at P,’ is 
given by 


Yi H gla ss Yi 
sin BP’N: = ————- sin N’P’P, . (27) 


Let the coordinates of a point B on the incident ray be xp, yp and 
Zp- The length of the segment BP.’ is given by 


——_ 1 


BP; = d, = [(xp - xp)’ + (yp - yp)’ + (zp - zp)*)2 , (28) 
x en 
and, COS @,=- -d d _- ’ ) 
3 | 
Yep ~- Yp , 
on 5, *§ = , (29) 
2 d, 
— = & | 


cos y, = q . 


3 


| 
_—— 


where cos a,, cos B, and cosy, are the direction cosines of the 
incident ray BP,’. 


—™. 
The angle of incidence at P,’ is BP,’N; which was obtained in equa- 
tion (27). 


, /s > > . + 2 
cos BP,’N; = cos @, COs an, 


(30) 


+ cos 8, cos By ,+ COS y, COS yn2. 


Replacing the direction cosines by their respective values defined in 
equation (29), 
d,cos BP’N} = (xp - Xp) cos o%,, 


(31 
+(yp - yp)cos By,+ (Zp - Zp) cos yp, 


which is a first function of xp, yp and zp. 
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The principal normal to surface S, at point P,’ meets the inside 
surface S, at point N,’. Repeating the procedure followed to obtain 
equation (7), the coordinates of N,’ will be 


, 


XN = Xp - t cos Are 
YN = yp - t cos By,, (32) 
and, Zn = Yp - t cos yn2- 
Points P,’,P, and N/ determine a plane which contains the incident 
ray BP’, and therefore point B. 
Again, taking the equation of the plane in the form 
A,x+B,y +C,z+1=0, (33) 


and substituting the coordinates of each known point, equation (33) 
gives 


A,Xp + B,yp +C,Zp =-1, 
A,Xp, + B,Yp, + C,Zp, = -1, (34) 
and, A, Xn + B,YN +C, Zn =-1. 


The unknowns A,, B,and C, can be determined, and 
A,Xp + B,yp +C,Z_p +1=0 (35) 


becomes the second function of Xp, VYpandZ,. 


The third function necessary to determine the coordinates of point B 
will be obtained by selecting an arbitrary vertical plane in the en- 
vironment where the deviation defined earlier will be measured. 
For convenience, this plane (P) is chosen perpendicular to the pro- 
jection, in the horizontal plane, of the observation line and at a dis- 
tance *d” from O’. Therefore O’M = d, and the coordinates of point 
M in plane (P) are 





Xx, = % + d.cosA, 
Yu = ¥ - d. sin A, (36) 


and Zz 
The direction of OM is given by 


cOS @,,/ = cos ' 
cos f,,! = cos (90°+ A) =-sinA, (37) 


and, COS y,," = cos 90°=0. 


~: 





m 





NON | 


Side 
tain 


dent 


36) 


— 


OPTICAL QUALITY OF WINDSHIELD DESIGNS 35 


The general equation of the vertical plane, in terms of its distance 
from O’ is 


(x - xX.) cos ayy’ + (y - y,) COS Bay +(z - Z,) COS yyyr- d=0. (38) 
The coordinates of point B will satisfy this equation and, 
(x _ X 9) cos A + (Yp- Y,) (- sin A) -~d=z9 (39) 


becomes the third function of xp and y,. 


Equations (31), (35) and (39) yield the values for xp, ypandzp. 


The intersection of the observation line with the arbitrary plane 
(P) at point B’ — xp’, yp’ and Zp’ is given by 


Xpr = Xx,» 
YB’ = YM> (40) 


and, Zp’ = Z,+d tanE. 
The deviation, defined as BB’ is therefore determined since the 
coordinates of points B and B’ are known. 
Two components of the deviation vector are of interest, 


AX =) ,- X 
B B’ 
ms (41) 
and AZ = Zp - Zp, } 


The component Ay cannot be used because it is a function of the 
selected arbitrary plane (P). The components Ax and Az can also 
be expressed as angles AA and AE respectively, or combined to 
give a magnitude and an angle. 


CONC LUSION 


To obtain a measure of the optical quality of a proposed wind- 
shield design, the ray-tracing development is repeated ata sufficient 
number of surface points for various locations of observation, The 
integrated effect of deviation vectors over the entire windshield area 
can be represented, for each observation point, by deviation patterns 
which can be used for analysis. 

The development given in this paper is in its simplest form. 
More elaborate techniques have been developed to provide for the 
study of binocular effects. 
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Fuel Injector Design Problem 


By Donold F. Zaowada 
Engineering Research and 
Advanced Product Study Office 
Ford Motor Company 


INTRODUCTION 


Designing a device which is described by differential equations 
with variable coefficients which are non-linear, non-analytic func- 
tions of one or more variables presents a difficult problem to a de- 
sign engineer. Although it is often time-consuming and costly to use 
a design procedure in which successive small changes are made in 
each and every parameter until the desired result is obtained, when 
one considers that non-linear equations are not readily solvable and 
that methods for synthesizing even linear systems are not too well 
developed, it is not surprising that a trial-and-error method is most 
often used for designing such a device. 

In some cases the number of trial solutions which need be made 
may be reduced considerably by assuming that the effects produced 
by one component part of the device can be isolated from the effects 
caused by the other parts, When an independent causal relationship 
can be established in this manner, at least one part of the device 
may be designed comparatively easily and the work of completing 
the rest of the design may be greatly simplified. 

This paper describes how this method of attack was used in 
designing a fuel injector. 


DESCRIPTION OF THE FUEL INJECTOR 


The principal parts of the fuel injector are shown in Figure 1. 
The problem was one of designing the cam and the ports so that fuel 
would be injected into the cylinder head of an internal combustion 
engine at a rate which was a predetermined function of the engine 
speed and time. 

At the start of each cycle of operation, the engine-driven cam 
forces the plunger downward. As ports A, B and C, which connect 
to the supply port, are closed by the plunger, the pressure of the 
fuel trapped in chamber C rises. Just as the last of these three 
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ports closes, the delivery port opens. The resulting surge of fuel 
from chamber C through the passageway in the plunger and int 
chamber M causes a rapid buildup of pressure in that chamber. The 
spring-loaded poppet valve is forced open and fuel flows through the | 
nozzle and into the combustion chamber. 

The downward movement of the plunger maintains the nozzle 
flow at the required rate until port B is opened to chambers C and M 
through the passageway in the plunger. The pressure in these 
chambers then falls rapidly, the poppet valve shuts, and flow through 
the nozzle ceases. Further rotation of the cam returns the plunger 
to its starting position and the cycle is repeated. j 

: Figure 2 shows the flow rate at which it was desired that fuel be 
injected through the nozzle. Because a constant quantity of fuel was 
to be injected during each cycle of operation, the flow rate required 
was proportional to engine speed and was normalized as shown in 
the figure. The curve has not been plotted beyond 125° of cam ro- 
tation because it was desired only that flow end as quickly as possi- 
ble at this point. The size of the nozzle and the rate of flow through } 
it were determined from factors relating to the combustion of the 
fuel which need not concern us here. 
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In designing an injector of this type to meet flow rate specifica- 
tions, two major obstacles are encountered, First, it is possible to 
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make the injector deliver fuel at the desired rate at only one engine 
speed. Second, the flow characteristic obtained at any engine speed 
is determined principally by the size and location of the delivery 
port and port C, and by the shape of the cam. These system param- 
eters are specified with difficulty because they do not perform their 
functions independently of one another. 

The first difficulty arises because, as the engine speed is in. 
creased, the increases in the pressures within the injector are pro- 
portionately greater than the change in speed. This causes a larger 
proportion of the fuel displaced by the plunger to be compressed in 
chambers B, C, and M, As a result, a smaller volume of the dis- 





| 


placed fuel can flow through the nozzle and the flow rate cannot vary | 


linearly with engine speed in the manner desired. 
To minimize this effect, the cross sectional area of the plunger, 
the approximate location of the delivery port, and the total plunger 


stroke were designed so that the volume of fuel under compression 


in chamber C would be as small as was practicable at the time the 
delivery port opened. Likewise, the volumes of chambers B and M 
were made as small as possible. These quantities were specified at 
the beginning of the design with no difficulty. The rest of the system 
was then designed to produce the desired flow at a particular engine 
speed chosen so that the deviation from the desired flow was within 
acceptable limits at other speeds. 

The difficulty encountered in designing the ports and cam results 
from the fact that since the delivery port cannot open instantane- 
ously, for all practical purposes, no cam can be devised which will 
cause the plunger to displace fuel at a rate sufficient to produce the 
sudden jump in nozzle flow required at 91° of cam position. This 
jump is approximated by greatly compressing the fuel in chamber C 
so that there is a sudden surge of fuel through the delivery port as 
it opens. Thus the flow through the nozzle is a function of the 
manner in which ports A, B and C close and the delivery port opens, 
and the shape of the cam after and immediately before the delivery 
port opens. 

The design problem is simplified by assuming that the fuel 
which flows through the nozzle is obtained from two independent 
sources: the expansion of the fuel compressed in chamber C before 
the delivery port opens, and the displacement of fuel by the plunger 
after the port opens. 

Considering only the time after the delivery port opens, for 4 
short interval the fuel compressed in chamber C expands into 
chamber M, At the same time, of course, most of the fuel displaced 
by the plunger also flows through the delivery port. After a time 
the fuel originally compressed in chamber C ceases to expand anda 
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state of dynamic equilibrium is reached in which the flow through 
the nozzle becomes dependent solely upon the rate at which the 
plunger dispiaces fuel. The time required to reach this equilibrium 
condition depends to a large degree upon the rate at which the de- 
livery port opens. If this port could open wide instantaneously, the 
equilibrium would be established virtually instantaneously. 

The effect of the initial buildup of the pressure in chamber C, 
designated hereafter as P., diminishes rapidly with time as the de- 
livery port opens. The equations describing the injector can be set 
up and solved under the assumption that this initial pressure buildup 
has no effect whatsoever on the nozzle flow. It may be expected that 
if a cam profile could be found which would yield a solution in which 
the desired nozzle flow was attained under such a condition, then 
this cam would also produce the required flow under normal oper- 
ating conditions if P. could be made to attain some value which 
would cause the flow due to the expansion of the fuel compressed in 
chamber C to make up for the deficiency in flow due to the plunger 
movement during the time the equilibrium condition was being 
established. 


DESIGN PROCEDURE 


Using the procedure outlined above, the profile of the cam was 
designed from the 91° to the 125° positions by assuming that the de- 
sired flow had been established and was being maintained solely by 
the movement of the plunger. 

To accomplish this it was assumed that the cam had rotated 
through 91°, that ports A, B and C were closed, that the delivery 
port was fully open, and that the pressures in the various chambers 
were initially those required to cause the desired flow through the 
nozzle at the start of the injection portion of the cycle, i.e., at 91° of 
cam rotation. Under such conditions, the pressure drop across the 
delivery port would be infinitesimal and chambers C and M could be 
regarded as one, 

Equation 1, the continuity equation, describes the system, but 
must be modified before being used. 


mg Sco oe rack atten 
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The symbols are defined as follows: 
A, - Area of nozzle orifice 


- Cross sectional area of plunger 


on 
1°] 


- Bulk modulus of fuel 


- Flow coefficient for all ports and orifices 


a 
2 


Pressure in chamber B 


7 
' 


Pressure in chamber C 


- Pressure in the combustion chamber (a constant) 


ie) 


- Pressure in chamber M 


= 
~ 


- Pressure drop across the nozzle 


be 


- Supply pressure (a constant) 


W 


vu yw wevweweanw > 
' 


a+ 


- Pressure drop across the poppet valve (a constant) 
t - Time 
V, - Volume of chamber B 


V..7- Volume of chamber C with plunger at top dead center 


Vim - Volume of chamber M 


X - Plunger position measured from top dead center 


As can be seen from the equation, it has been assumed that fric- 
tion in passages and orifices, inertias, head losses, transmission 
delays, acoustic effects and leakage are all negligible, and that the 
nozzle is a sharp edged orifice. It can also be assumed that the 
supply pressure and the combustion chamber pressure are constant. 

The poppet valve is designed sothat when it is open the pressure 
drop across it will be virtually constant, regardless of the variation 
in the flow through it. Since the pressure in chamber M rises 
rapidly, no great error is introduced by assuming that this valve 
opens instantaneously when the differential pressure across it 
equals the constant pressure drop across it when it is open. Thus 
Pi- Pp < P when the valve is closed; P,,, - PR, = Bh when it is 
open. 

When the valve is closed, P}, equals P,. Asa result 
(3A) P, =0 

dP, when P< P, + P 


(4A) i = 0 


t 














.DA 
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When the valve is open, P;, equals P,,, - P,;. Thus 
(3B) P,, =Ph- ht - P, 
dPp dP, when P,, > P, + P; 


If itis assumed that fuel is flowing through the nozzle, the poppet 
valve must necessarily be open and equations 1, 2, 3B and 4B can 
be combined to form equation 5. 

aX V. "Se A,X + Vm bai Vb dP, 
dt B dt 








) + Ay Cy /Pm- R-P, 


The last term on the right side of equation 5 is the expression 
for the rate of flow through the nozzle. This expression was equated 
to the flow rates desired at a particular engine speed, and P,,, was 
computed and plotted as a function of time. The first derivative of 
P,, was then used as a forcing function in solving equation 5 for X 
and its first derivative on an analog computer. Having obtained X in 
this manner, the shape of the cam profile from 91° to 125° of cam 
position was known, 

Because a portion of the fuel displaced by the plunger remains 
in compression and never reaches the combustion chamber, as was 
noted before, a different cam profile was obtained at each engine 
speed at which equation 5 was solved. Figure 3 shows plunger ve- 
locity curves obtained at four different engine speeds. The velocity 
of the plunger has been normalized by dividing it by the velocity of 
the cam. 

To complete the design a plunger velocity curve was arbitrarily 
plotted from 125° to 360° of cam position. Since the only important 
functions the cam performs inthis portion of the cycle are to reopen 
port B as a means of ending flow through the nozzle, and to return 
the plunger to its starting position, the exact shape of the cam in 
this region was not critical. Likewise, the velocity curve from 0° to 
91° was plotted arbitrarily, care being taken that the plunger would 
be moved down far enough at 91° so that the volume of fuel in 
chamber C would be equal to that specified at the beginning of the 
design. The complete cam profile was then obtained by integrating 
the plunger velocity function. 

From among the family of cam velocity curves which can be 
plotted, one must be selected which will produce the optimum re- 
sults over the entire range of speeds at which the engine is to be 
operated. Ordinarily it is desired that the percentage deviations 
from the desired flow characteristic be approximately equal in 
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Figure 3. Plunger Velocity Curves Obtained by Solv- 


ing Equation 5 at Four Engine Speed 


magnitude at the maximum and minimum engine speeds. Following 
a rule of thumb, this result can be obtained by solving equation 5 at 
the engine speed which is the geometric mean of the maximum and 
minimum speeds. When some other result is desired, several de- 
signs of the injector must be completed and evaluated. This was 
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done in the design described here and the cam obtained at 3000 RPM 
was judged to perform best. 

It was necessary to return to the normal conditions of operation 
to determine what amount of pressure buildup in P. would cause the 
nozzle flow to be the desired value during the time the dynamic 
equilibrium condition was being established. 

The equations for the complete, normal cycle of operation were 
written by equating the rates of flow into and out of chamber C, and 
into and out of chambers B and M. The assumptions regarding the 
poppet valve were used again to simplify the solution. 


Veo - ApX (dP 
i ie ) + Ae Cy J/P.- Ps +AgCy /Pco - Pm 


(6) AL at B dt 
(TA) AJC, /P. - = 3 Ca when P,, < P, + P, 

—_— — V. bd dP. —_ =m — 
(7B) AgCg/P. - Pm =— +AnCg /Pm - Pt- Py 


when am > PR + _- 


The symbols are those used previously, with the following ad- 
ditions: 


A. - Area of delivery port (a function of X) 


A, - Combined area of ports A, B and C, and the effective 
area of leakage paths in parallel with ports A, B and C 
(a function of X) 


It was assumed that the portsfunctioned as sharp-edged orifices. 
Since leakage to the supply port was small and could affect the so- 
lution only during the short time all of the ports were simultane- 
“usly closed, or almost closed, no great error resulted from repre- 
senting the leakage paths by a sharp-edged orifice of constant area. 

The plunger velocity function obtained at 3000 RPM was used as 
a forcing function in solving equations 6, 7A and 7B on an analog 
computer. The pressure buildup required to produce the desired 
flow rate was most easily determined by trial and error variation of 
the locations of the ports. This was not a very lengthy process. 
Ports A and B were designed to close before port C began closing. 
These two ports and the delivery port were also designed to be 
fairly large in area; the latter, in order not to impede the flow of 
fuel at the start of the injection, A and B, in order that chamber C 
could refull rapidly (through paths not shown in Figure 1) during the 
return stroke of the plunger. The area of port C was arbitrarily 
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specified to be a value which could not produce a large pressure 
buildup before port C began closing. 

Thus the problem merely became one of varying the locations of 
port C and the delivery port until a flow rate approaching the de- 
sired rate was obtained at 3000 RPM. A computer record of the 
result is shownin Figure 4. Figure 5 shows flow curves obtained 
from the computer at four different engine speeds. In spite of the 
many simplifying assumptions which were made in constructing the 
mathematical model of the injector, data obtained from laboratory 
tests of an injector correlated quite well with data obtained from 
the computer, 

The assumption, without prior justification, that the effects pro- 
duced by some parts of the injector could be isolated from the 
effects produced by the other parts greatly reduced the time re- 
quired to complete the design. 
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Analysis of Spring Installations for Low-Frequency 
Suspension of Vehicles 


By Raymond R. McHenry 
Engineering Staff 

Ford Motor Company 
Dearborn, Michigan 


INTRODUCTION 


The primary ride frequencies (bounce and pitch) of modern U.S. 
automobiles are in the range of 60 to 75 cycles/minute, and the 
damping level for these oscillations is approximately 25 to 30 per- 
cent of critical. Moderate reductions in frequencies (to the range of 
40 to 60 cycles/minute) without critical damping can, due to the 
slower resonant oscillations, produce a ride of a floating nature 
near the sea-sickness range, in addition to creating serious prob- 
lems with vehicle height and attitude for static load changes and 
dynamic weight shift. If load-leveling and some form of compen- 
sation for dynamic weight shift are incorporated, the primary ride 
frequencies should logically be reduced to the range where the 
damping-force level of the conventional “shock absorber” combina- 
tion of blow-off valves and hydraulic orifices will provide critical 
damping (approximately 30 to 40 cycles/minute). A modified form 
of damping-force diagram has recently been developed for low- 
frequency suspensions to more readily provide critical damping of 
ride motions and to eliminate the strike-through problem at the 
ends of wheel travel (the damping-force is a function of both velocity 
and position of the wheel relative tc the vehicle), 

In a current-production vehicle the change in wheel load per 
inch of wheel travel is approximately 10 to 12 percent (referred to 
sprung portion of wheel load), For a low-frequency vehicle (30 to 40 
cycles/minute), the corresponding change must be reduced to ap- 
proximately 3 percent per inch of wheel travel. For this reason, the 
geometry of the spring installation is much more critical. 
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I. CALCULATION OF WHEEL RATES INCLUDING 
THE EFFECTS OF SPRING INSTALLATION GEOMETRY 


Definitions: 


L._. = Wheel load (lbs.). 
L. = Spring load (lbs.). 


x = Change in spring length from design length (inches, posi- 
tive in direction which increases spring load), 


y = Vertical travel of the tire-to-ground contact point, relative 
to the vehicle, from design position (inches, positive for 
movement toward vehicle), 


Neglecting dissipative forces, the work done in moving the wheel 
from its equilibrium position is equal to the potential energy stored 
in the spring by the corresponding change in spring length, Ex- 
pressing this in differential form: 


L,, dy = L, dx (1) 
The rate of change of wheel load with respect to wheel travel 
dL, 
(wheel rate) is defined: Wheel Rate = dy 
From equation (1), 
dL a( =) . d x =(2) » 
——w = —(L_—) = L + 2 
dy Say) © “s ay* “4 


The spring rate along its own center line must be de- 


dL. 
, ax , 
termined from the characteristics of the spring medium. 

x = f(y) in some cases can be directly derived from an analysis 
of the spring-installation geometry (e.g., simple trailing or leading 
arm suspension). If a functional relationship cannot be derived 
readily, corresponding values of x and y can be obtained from a 
digital computer program of the suspension and spring-instaLlation 
geometry and an empirical polynomial equation can be fitted to this 
data. In either case, when a functional relationship is obtained be- 
tween x and y, the first and second derivatives of x with respect 
to y can be calculated for use in the expression (2) for wheel rate. 

The first term of expression (2) is the *geometric rate” due toa 


dx 
non-linear relationship between x and y. A linear relationship (= 


d x 
= constant, of. 0) would eliminate the first term and reduce 


dy 
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expression (2) to the generally used approximation for wheel rate 
: ; , dx 
(spring rate times nominal installation ratio, - squared), The 


first term becomes extremely important for low-frequency sus- 
pensions since the required wheel rates produce only a 3 percent 
change (approximately) in wheel load per inch of wheel travel (e.g., 
1000 lb. wheel load, 30 lbs./inch wheel rate). In such an installation, 
the assumption of a linear relationship between x and y in the cal- 
culation of wheel rate for a nominal 0.500 ratio with arate of change 
of ratio of 1.5 ‘%/inch (+ .0075 inch) at design height would produce 
a + 50 percent error in the effective design height wheel rate: 
From equation (2), 


the design height wheel rate = (2000)(7 .0075) + (120)(.250) 
30 7 15 lbs./inch 


Nominal + 50% 


Since detection of a 1.5 percent change in ratio (dx/dy) per inch 
of wheel travel is almost impossible by layout drafting methods, 
low-frequency suspension installations must be analyzed by mathe- 
matical methods. The complexity of the functional relationship be- 
tween x and y will in most cases require a digital computer pro- 
gram to make the analysis practical. 

The use of air bags for a spring medium complicates the analy- 
sis by introducing a variable load-supporting area. For this reason, 
the analysis for a gaseous spring medium will be demonstrated with 
a hydro-pneumatic strut. 


Il. CALCULATION OF WHEEL RATES 
FOR HYDRO-PNEUMATIC SPRING MEDIUM 
Assumptions: 


1) Ideal gas laws. 


( 

(2) Isothermal leveling process. 

(3) Polytropic process for dynamic spring action. 
( 


4) Precharge pressure measured with gas at ambient tempera- 
ture. 


P, V, 
Static Gas Volume, V, = —‘p,~ in? (3) 
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Where P, = Precharge pressure (psia) 
P = Pressure under static load (psia) 
Vv, = Precharge volume (in.° ) 


From the polytropic process for dynamic spring action: 








© G- psia (4) 
V, = (V, - Ax) in” where x = strut travel from design 
length (inches, positive 
in compression). 
A = load supporting area (in.? ) 
: k = polytropic exponent 
V k P, 
= — - "oe 3ic 5 
: m8 (Gres) Ze ° 
~ V 
1 
PA 
Strut Load, L, = (P, - 14.7) A = 7 waak ~ 147A Ibs. (6) 
i 4 
(1-7) 
1 
dL; k Py A’ 
Strut Rate == Ib./in, (7) 
eeu SUA, Ax\k+1 
*70-8) 


By substituting equations (6) and (7) into equation (2) we obtain: 








dL, P,A d*x , 
heel Rate, ——# = | ——1+-"—- - 14.7 : 
Wheel Rate, —7, | , . ARy 4 | l ay? 
Vi (8) 
- aan | (ey lb. /in 
7% Ax \¥ dy / 
LY, ( -¥F ) 


In conventional independent front suspension spring-installations, 


d“x , ; : 
the 1? term in equations (2) and (8) is negative throughout the 
C 
. strut travel ; 
range of wheel travel (ratio of —————— larger in extended po- 
wheel travel 





sition than in compressed position). In conventional trailing-arm 
rear suspensions, the term is positive throughout the range of 


ANALY! 


cravel. 
observe 


(1) 


(2) 





4) 
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travel. The following effects are indicated and have actually been 
observed on experimental low-frequency vehicles: 


(1) As the primary ride frequencies (bounce and pitch) are pro- 
gressively reduced by decreasing the spring rates, the front 
suspension becomes unstable. 


This instability is due to a non-uniform geometric 
effect. The *geometric rate” is negative throughout the 
range of travel; however, it is larger near the extended 
(rebound) end of travel and it therefore produces a nega- 
tive effective wheel rate in this region of travel when the 
spring rates are sufficiently reduced. 


(2) With gaseous spring media, the variation between static and 
dynamic spring rates (due to change of the polytropic gas 
exponent) produces an increased percentage variation in 
effective front wheel rates and a decreased percentage vari- 
ation in effective rear wheel rates. 


A change in front spring rate appears as a larger per- 
centage change in wheel rate since the effective front 
wheel rate (equation 8) is the difference between the term 
involving spring rate and the “geometric rate.” 

The opposite is true on a trailing-arm rear since the 
effective rear wheel rate is the sum of the term involving 
spring rate and the *geometric rate.” 


It is readily seen from the described effects that positive values 
d*x : ; , 
of dy? are desirable from the point of view of exponent changes; 


however, gas volume requirements dictate that the positive values 
must be small, 
Advantageous application of geometric effects (large negative 


d x ' ‘ 
values of dy? to reduce wheel rates) requires a spring medium with 





a rate which is insensitive to excitation frequency. For the follow- 
ing, a conventional round wire extension or compression spring will 
be considered. 


Il. REDUCTION OF EFFECTIVE WHEEL RATES 
BY MEANS OF SPRING INSTALLATION GEOMETRY 


(1) Geometry required for linear wheel rate from linear spring: 





In equation (1), let Lw = Lwo + Kwy, Ls = Lso + Ks x and 
integrate: 

























* #86 #2 
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2 2 
x Ky 
- + L sx wy + Lwoy (9 
Where 
K. = Spring rate (lbs./in.) 
K. = Wheel rate (lbs./in.) 
a = Static spring load (lbs.) 
L = Static wheel load (lbs.) 


The locus of this equation is a hyperbola with center a 


Lso Lwo 
*-= 3° - K or the degenerate case of a pair of intersectin 


Ww 





Lwo 2 
lines if Ky = K;\7 ) : 


Solving equation (9) for x and differentiating with respec 


to y: 
ae ? wa) L \: 
a + 2 , me 
K. N(R “) y g K, ii a " 




















x — 
ox _ Lyo + Ky 
dy L,, + K,x 7 
dx 2 
d*x _ K,,- K,( dy (12 


dy? ss Lig) + KX 


(2) Requirements for load-leveling: 





2 
If large negative values of ra are utilized to reduce effec- 


tive wheel rates, the rates of a system with load-levelix 
and a fixed relationship between spring and wheel travel will 
be sensitive to loading. The rate reduction effect will in- 
crease with the static spring load (see equation 2), and sym- 
metry of rate-deflection will be altered by load changes. 


If critical damping can be obtained for the driver-only con- 
dition, the optimum rate vs. static-load relationship woul 
appear to be a compromise between maintaining critical 
damping (rates decreased as load is increased) and mini- 
mizing the effect of friction on height stability. Such a com- 
promise will require determination of the actual form and 
level of damping as well as measurement of friction. I 


ANALY 





Spring 


Sprin 
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order to demonstrate the relationships involved in a system 
with load-leveling, it will be assumed for simplicity that 
constant ride frequencies are desired in the following calcu- 


lations: 
fg \? , 
*.. * Lo (Gans Where f, = Design height 
frequency (13) 
(cycles/ minute) 


By substituting equation (13) into equation (12) and letting 
x = 0, we obtain: 


CG), ESL. GHs)- G)] os 


The required relationships of the indicated parameters are 
shown graphically in Figure 1. 


Utilization of geometric effects for reduction of effective wheel 


rates is advantageous from the viewpoints of efficient application of 
the energy-storing medium and package considerations. 


IV. REDUCTION OF REQUIRED VOLUME OF SPRING STEEL 
BY MEANS OF SPRING INSTALLATION GEOMETRY 


Spring Energy at Design Position: 











P = Spring load (lbs.) 
i vt P Ks = Spring rate (lbs./in.) 
> "ae . Sagith 
Lp = Design wheel load (lbs.) 
Pp , “a ; (sering travel ) 
Ks Rp = Design ratio me ee 
1P° Lp 
S > 2 2B asain 
pring energy = > Ke’ Rp 
1 LD\? 
Design position spring energy = — (15) 
2Ks 








4 Energy in Spring = A Energy at Wheel: 





By means of ratio variation effects, the effective wheel rate pro- 
duced by a linear spring can assume a variety of forms. For this 
study, two forms of wheel rate will be considered: 














| Ks = 500 Ibs./In 


DTT #4 = 36 cpm 









+20 saa 
(X-Xo)o | 
(Inches)-+ 












let 
E 
r 
— ~--- - let 
80 
ria n- Mz 
(1) K, = K,,. = Constant (produced by a hyperbolic relationship 
between spring and wheel displacement including 
the degenerate case of a hyperbola which is 4 E 
n 


straight line). : 
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(2) Kw =Kywot+ Ky y* = Symmetrical rate build-up (parabolic 
rate curve). 


The two cases will be covered by solving for case (2) and 
letting K,,, = 0 to produce case (1). 


K, = K,, + K,, y’, for 4.5 inch maximum wheel travel: 
_ = K a 20.25 K. 


Kum (= ) Kwm/ Kwo™! 
K., = 1 + 20.25 K. ’ let 20.25 =a, K. 1= a K..., 


) 





) 


Kk, = K,,{i + ay") (16) 


Wheel load = Lp + K,,,, | “(l+a y”) dy 





— 
A Energy = | (Lp + Kwoy +—s—Yy’) dy 
it ’ (18) 
Ko Bue @ y* 
= Lpy + 7) y 4 “— 
1 (DY KwoY 
Maximum Spring Energy = 2K Ry +Lpy +9 
: ‘ (19) 
Ky, @ Yy 
ra" 
Ewe L, 
let B = ———, 8 ®* =—~ 
. K.,(R,)? K 
E = 52 ( 36 2y + y_ + e y-) (20) 
max D) PV - A 6 5 
let y_.., = 4-5 inches 
Lp 20.25 a 
Maximum Spring Energy, E., ~. 3 Bb +9 +- —— 68.344-—) 
Lp 187.5 ? 
§=—— ———) 
Kwo f , 
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Active Volume Required in Round Wire Compression 


and Extension Springs: 





























2 K r =z 1 p* 
oh 
f= Pry 
— 
K 25° som ‘ a ; 
f 8ND got B64) HCz) --se (22 
p = 78d 2K ‘2’\64D" “Ts 
8D 
- < a? A2 aT 
Active Volume = (N) (7 D) (=2\ = — 
2 4 
Energy _ inch-lbs., V (in.’) = (e E (23) 
r 2 
Active Volume Required = “cre[s + B (482.8 
a D 


+ (5.76 x 107*) (f,2) (1 + 3.38 a)| 


Let G = 11.5 x 10°, Lp = 1,100 lbs., S (uncorrected) = 110,000 psi 
87.5\" a ; 

V = 2.09 [9 + B (415) + (5.76 x 107*) (f, ) (1 + 3.38 a)| in.* (24) 
D ; 


Where f.. = Design height frequency (cpm) 


K y ) 

—. 
eK (R,)? 

, =-1l 
a K A 

ee a 
20.25 

Figure 2 shows the effect of variation of 8 on active spring steel 


volume requirements for the following two cases: 


(1) Constant wheel rate (a = 0). 


(2) Parabolic wheel rate with rate at full travel double that at 
design height (a = 0.0494), 





AN 
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Design and Calculation of Progressive 
Compression Springs 
By Albert E. Roller 


Pontioc Motor Division 
Generol Motors Corporation 


1. INTRODUCTION 


This analysis deals with round wire helical compression springs 
in which a progressive rate is produced by means of a varying 
number of active coils. With the derivation of a simple equation by 
integration it is possible to calculate the rate of such a spring. 






Vtarting of progrevtive rale 
A 


nga 


_ 


t 


deflection f- 





In most ordinary springs, deflection is directly proportional to load, 
Sometimes, however, progressive rate is a requirement. It is pos- 
sible to design a compression spring with a definitely increasing 
rate (e.g. by using a variable pitch of coils). During deflection of 


6! 
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such a spring, touching of coils starts at any point depending on the 
spring design. This reduces the number of active coils and in- 
creases the rate of the spring. A spring working under these con- 
ditions may be called a progressive spring. 


Symbols. 


G = Torsional modulus of rigidity lb.in~ 
d = Diameter of wire inch. 
D = Mean diameter of coils inch, 
n = Full number of active coils before touching of 
coils 
n, = Number of eliminated coils at deflection 6,, 
n, = Number of active coils depending on deflection - 
c = Rate of spring lb.in 
§ = Deflection of spring inch. 
5» = Deflection of spring between progressive limits inch, 
P = Load on spring lbs. _, 
S. = Shear stress lb.in 
k = Stress factor for spring design depending on 
=. 
ratio — 
d 
tT = kx S&S, = Resultant shearing stress when stress- - 
factor applied lb.in 
. _ 304°. : 
f = —/— = Frequency in cycles per min when ap- 
plied load on spring is P=mxg 
g = Acceleration due to gravity = (12)(32.2) inch,sec” 
s') sascha 
. “a S* Factor for simplification 
3 
b = 8D x 2p = Factor for simplification 
dG 5p 
W = Work done on the spring inch.lbs. 


Feature of progressive springs. 
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If the frequency f = a’ = — is to be constant in the range of 


 @ 
the progressive spring rate, the ratio = must be constant. 


From this results the derivation, 


2.2.6 
m d6 P 
_g ;@P_g P 
| o= 5 f pag p 
P=-P xe ® (1) 


Equation (1) indicates that the spring rate follows the law of the 
exponential function 


y = eax 


In this function all curves intersect the ordinate y at unit 1, 





0 He 7 * en 
S)y 


igure 3 (for — = const. only) 


2. General Relations for a Progressive Compressive 


Round Wire Helical Spring. 


From the known relation for this type of spring, 


3 
a5 = He sn, aP 
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: —_ dP ' ; 
it will be recognized that a deviation of a5” the rate of spring, is 
) 


only possible by changing the number of active coils. This can be 
arranged by producing a spring with variable pitch of coils operating 
as explained in the introduction. The following analysis provides a 
method of computing progressive springs. 


If the design of a spring provides for the elimination of active coils 
in proportion to deflection 6, the number of eliminated coils will be, 


§-— (See Fig. 4) 


NMa=n- 6— (2a) 


Point A in fig. 4 where coils begin to roll off may be determined 
according to the requirements of the spring application. 


Note, Considering fig. 4 the spring calculation has to be made in 


parts. For range O-A the ordinary spring equations are ap- 
plied. For range A-B the following equations can be used. 


Development of the Basic Equation. 





Substituting the value of n, from (2a) into (2) we obtain, 


8 D* | n,] 8 D° 8 pD* nn...) 
dé = — n- 6 P\dp =|—— n - —P §idP 
d*G | 5p | \d*G = = dG by 

a b 
dé = (a - bé)dP 
P= | bs. + mis - bs) + C 
in. - ilies 


For point A we get, P = O and 6 = O therefore 


O see iar 
b 
ae 
C =, ima 
p =+{ma- in(a- bs)] = + in _ 
» * b 1 dD, 





Fr 


Fu 
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cia ' b 
P =; [In 1 - In(l - = 6)] 
— b 
pPe- 5 BG -5 9) (3) 


This is the basic equation for progressive spring rate between 
point A and B in fig. 4 where, 





b n 8 D® on 
o- s — ~ ¢ b =- a. 
a = §°SS ss 


Conversion of the basic equation shows that it is possible to solve 
mathematically, 














I 
e-bP zi. ? ray 
a 
[ a rs 
og Vee (3a) 
a J 
b 1 
Forn=n it will be - = — 
P a 6p 
: , b n, : 
For 6 = dp itwillbe — x 6 = —® and we find 
1 ] 
P=--— ln ( - Rp) (3b) 
b n 
From (3b) we get, 
e~bP = 1 a Np 
n 
Np a 
ae re 
= 1 - © - 
Furthermore, 
a, 
- bP x loge = log(1 - P) 
8 D* n n 
-—,— --—2. P- loge = log(1 - —2) 
d‘G Dp 5 n 
8 log | 
| 5p = --He . Np - Pp a e log e = 0.43429 (4) 


ns | 
y —— | 
log(1 - ) 


- a 
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Note, For elimination of active coils in proportion to deflection 6, 
"any other function can be used. In the case that a solution by 
integration is impossible a graphic method usually will be 
satisfactory. 


3. Stress Formula, 
Computing of stresses, of course, will be the same as in ordinary 
compressive springs. 
Therefore shearing stress at inside of coil, 
8 D 


els hee ee lb.in™* (5) 


k = *Wahl correction factor” 


Note, Within the progressive range do not calculate the stresses in 

~- proportion to deflection as often is done when computing 
ordinary springs. In this case the stresses must be calcu- 
lated by using the force P because P is not proportional to 
deflection, 


4. Work done on the spring. 


Considering fig. 4 again, we see that the work must be computed by 
parts. 


For the range O-A (triangle OAD) we get, 
WwW, = — om inch.1lb, where Pa = Force at point A, 


For the range A-B, 


Area ABC, 


5=5, 5, 
- 1 b 
W, = | Pdd' = | - = In(1 - — 6)d6 
- b a 
=0 0 
6; 
1 +f b 
Bee ln * 
b } (1 a 5)d6 
For 1 - - 6=V we get 





PR¢ 





nf 
e 
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On 

a ze , ae ee -_ Yp 

W, =p J mVv-dv = [v-mv- Vv}, 
a, b b b_,,5 
= om 4 aoe Sib « & @ « 18 « @ Bit 
pi = 5] n ( = 5) (1 5 O)) 








= 1 , 
2 a - D . =. . . mn ’ 
| cs {1 —. [In (1 ; 5... ) 1 1 (6) 
Area ACED, 
W, = Pa * 5, 


Total work done on spring therefore, 
W = W, +W, +W, inch. lb, 


Example, Design of a helical progressive compression spring with 
the following assumed data, 





Wire diameter d = 5/16” 

Mean diameter of coils D.. = 3” 

Number of active coils at point B n, = Np = 4,5 
Maximum load on spring at point B Pp = 500 lbs. 
Load on spring at point A P, = 200 lbs. 
Spring index rt = 6,4 

Stress factor k = 1.23 

Torsional modulus of rigidity for steel G = 11.5x 10° Ibin™* 


Solution, (See also Fig. 4) 
Maximum shearing stress at inside of coil, 


8 D 8.2 
max ae bie the , 


“ 


F-91053 * 500 x 1.23 = 102 600 Ib.in™ 


Spring rate at point B, 











4 4 6 
_d@P_ qd-G_ _ 0.3125 x11.5x10 _ wi 
Cp = a 8. D>. np 8x2? x 45 381 Ib.in 
Frequency 
30 1/c’_ 30}/381 x 12 x 32.2 _ ol 
st 1 Vs ; =h 500 . ao — 
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If equal frequency at point A and B is desired we obtain the spring Ps 
rate at range O-A, 2 
>, = Ao 966 See © 152.5 Ib.in.” 
C, =Cpz Sy 500 - -in. | 
Full number of active coils before touching of coils, 3 
™~ 

ar 2 CB 381 | — N 
n = gp; C. hy C, = 4.5 152.5 11.25 coils ‘ 
Therefore the number of eliminated coils at deflection 5, from 
point A to point B, 
n, =n-n, = 11.25 - 4.5 = 6.75 coils 
Deflection from O to A, 

Pa 200 — 

*. ms =" 
Deflection from A to B by using equation (4) . 
gD __loge __ 
bn = - aig n,(Pp- P,) 

log(1 - 2) 
8 x 2° 0.43 429 
es - Wr 29 in, 

Sp = "0.3125"x 11.5 x 10° * 7 x (500 - 200)F eons _y ~ 1-29 in , 
, a eee Fe 6.75 _ 1 | 3 
> “GG, 0.3125" x 11.5 x 10° 1.29 327 
b Np . 6.75 i 
a n- 6 11.25 x 1.29 2.15 


p 


Calculation of curve A-B, 


0 lbs. 
0.2 | 0.093 | 0.907 | 6.81014 - 6.90776 = - 0.0976 r 32 lbs. 


0.4 {0. 186 | 0 0.814 [6.70196 - 6.90776 =-0, et. 67 lbs. 


0.6 | 0.279| 0. 721 |6. 58064 - 6.90776 = - 0.3271 107 Ibs. 
0.8 10.372 | 0.628 [6. 44254 - 6.90776 = - 0.4652 | 152 lbs. 
=a 0.465 | 0 0.535 {6 .28227 - 6.90776 = =-0,6255| 205 lbs. 

| 0. 558 | 0.442 = .09131 - 6.90776 = - 0.8165 267 Ibs. 
1. aati 0.600 | 0.400 | 5.99146 - 6.90776 = - - 0.9163 | 300 Ibs. =| 
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Figure 4 
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' 


Work done on the spring see fig. 4, 





Regarding area OAD, 


Pa *5 200 x 1.31 





a -——« = = 131 in.lb. 


Regarding area ABC, 


_a 1 5, b b. 1 
WM =H pill- 7 Splint - 5 5,) - 1) +1) 
W, = 2,15 x 327 { [0.400] [- 0.9163 - 1] + 1 } 


W, = 165 in.Ib. 

Regarding area ACED, 

W, = Pa 5p = 200 + 1,29 = 258 in.lb. 

Total work done on spring during deflection O-B, 

Ww = W, + W, + W, = 131 + 165 + 258 = 554 in.lb, 


This example is based on the assumption that the frequency for 
points A and B is the same, 


Based on the equation outlined above most requirements in pro- 
gressive compression spring design can be calculated as practised 
in cOmmon spring computation, 
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Experiences with Highly Loaded Gearing 


By Forest McFarland 
Buick Division 
Generol Motors Corporation 


INTRODUCTION 


As a result of seeing its application for many years, I am very 
conscious of the importance of mathematics to engineering design. 
I hope in this discussion to relate some experiences in actual 
designs and developments which may be of value to you. 


GEARING IN LAY SHAFT AUTOMOBILE TRANSMISSIONS 


In the early thirties, attention began to be paid to making the 
gearing in automobile transmissions not only durable but quiet. To 
do this, engineers began to employ helical gearing. It was not 
enough simply to cut, harden, and assemble the gears — several 
rules had to be followed. I will describe the requirements for 
satisfactory operation, found as a result of the development of the 
first Packard helical gear transmission, in which I had a part, under 
the direction of W. R. Griswold, then Chief Research Engineer. A 
later version of this transmission employing the same principle of 
gearing is shown in Figure 1. 


1. The gears must be cut and “shaved” very accurately for 
lead, profile, and index. These particular gears were also lapped by 
the “Incolap” method. 

2. The profile of the gears must be accurately held, but pri- 
marily must be designed to provide “soft” contact at the start of 
engagement, which means modifying the involute. Figure 2 shows a 
typical profile deviation chart. You will note from the modification 
charts the action is from the pitch point on, or recess action. I have 
been told that the old clocks with wooden teeth employed mainly 
recess action. I have not seen it, but have mentally theorized on 
benefits of this action and suspect a good case can be made for it 
with gearing of questionable efficiency. 
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Figure 3. Rigid Ball Bearing Gear 


Mounting. 


3. The gears should have a helical contact ratio of two or more 
if possible. Helical contact ratio appears “worth more” than invo- 
lute contact ratio to reduce deflection and reduce *PV” values, 
Angles on the constant mesh gears of over 45° were called for in 
the forward constant mesh gears and angles of the same lead, of 
course, in the first and second speed gears, to balance the counter- 
shaft end thrust. 

4. The gears must be rigidly and accurately mounted. High 
helix angles make this problem more difficult. Figure 3 shows ball 
bearing mountings in two of the gears. This was a development all 
in itself. The completed transmission is again shown in Figure 1. 
This unit was extremely quiet and very rugged. It was practically 
impossible to determine whether the car was in second gear or 
direct drive except by engine speed. High helix angles make for 
more expensive mountings, but the additional cost was felt to be 
justified in this case. 


Since transmission gears operate at full engine torque only a 
small percent of their life, the question arose as to what stresses 
could be allowed, The most realistic work I know of was done by 
J. O, Almen and John Straub, published in the middle Thirties, when 
they were with General Motors Research. Their method, briefly, 
was to examine the service life of various transmissions, calculate 
the stresses by a modified Lewis formula and then run these units 
at full torque in first and second gear. The results showed that a 
transmission good for 100,000 cycles in first gear and 300,000 in 
second would give satisfactory life in service. 
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Figure 4, Rolls-Royce Aircraft Gears 


REAR AXLE GEARING 


Hypoid rear axle gearing, introduced by Packard in 1924, has 
become standard in the industry. In the design of rear axle gearing, 
we are again “plagued” by deflections which cause the contact to 
move particularly lengthwise on the gear teeth. As is true of many 
other parts, rear axle gear-set components are designed primarily 
for stiffness rather than for strength. The gears cut by the Gleason 
process provide deflection of the parts by a “crowning” effect of the 
teeth to keep the contact “on the teeth” rather than over the edge of 
the teeth. This applies particularly lengthwise on the teeth but also 
to the tip of the teeth. The more the deflection, the more must be 
the “crowning.” It is, therefore, highly important that the mountings 
be kept very rigid to keep deflection to a minimum; to minimize 
crowning; to keep Hertz stresses to a minimum; to insure capacity 
and also to keep spiral overlap to a maximum for quietness. De- 
flection tests are run on the gearing in drive and coast, and the 
results very carefully studied. Balance, of course, must be obtained 
in the design, It is desirable that the contact under no load be in the 
vicinity of the middle of the teeth and stay there under load. As a 
result of a careful analysis on the writer’s part of the effect of each 
individual deflection on movement of the bearing along the gear face 
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SLIDING OCCURS PERPENDICULAR TO LINE 
OF ACTION AT THE POINT OF CONTACT. 
THE VELOCITY VECTORS Vv, AND V, HAVE 
INSTANTANEOUS TANGENTIAL VELOCITIES 
AROUND THEIR RESPECTIVE CENTERS OF 
ROTATION AS INDICATED BELOW 

RPM(ORIVER) 21 p .y FT 

60 “et "se 

RPM(ORIVEN) 27 2 _y 

—e" o nV x 2°‘%2 
SLIDING VELOCITY = (VV) « 








035xNOR FORCE(|A+ h,) es 
FACE WIOTH (= +e.) 


fad 


DRIVER SURFACE STRESS: 





PV.=Px(VeV,)* Be. 


Figure 5. Surface Stress & PV Value 


of the particular gear set studied, it was concluded that stiffening 
the pinion shank, particularly at the rearward end, provided more 
control of contact for the money expended than any other one item 
of design. As a result, the largest rear bearing bore commercially 
feasible was worked into this design. This was provided even to the 
extent of obtaining a cone of special inside diameter. The results 
finally obtained from this design were gratifying, as the gear set 
was definitely smaller than several others, even though behind the 
largest engine in the industry at that time. 


AIRCRAFT ENGINE GEARS 


During the last war, the writer had the opportunity to observe 
the behavior of the gears (approximately eighty in number) of the 
Rolls-Royce Merlin Engine as produced by Packard. The most 
challenging problem was the elimination of “scuffing” —a combina- 
tion of pitting and scoring in the reduction gears which were very 
highly loaded, being approximately 4,000 lb.-per-inch of face at 
the start of production and ending up at approximately 6,300 lb.- 
per-inch of face on the final *dash 9” production engines, operating 
in an engine oil with no particular anti-scuff additives. 
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Figure ¢ . PV Value for One Set f Teeth in Action 


A tremendous amount of work was done on this gear-set. It was 
found that trouble-free operation required:- 


Proper profile of teeth, 

Accurate mounting of gears in good alignment. 
Parallelism of teeth within close limits. 
Proper surface finish on teeth. 

Adequate lubrication at right times. 


ue. © Ne 


Solution of this problem was very elusive as we could repeatedly 
prove that violation of any one of these items would not cause 
trouble if the other items were within proper limits. This was a 
case where pure “logic” would lead one astray. Taking sufficient 
liberties with any two or three items gave trouble. I am mentioning 
this point in particular as I have seen the same situation arise in 
other troubles resulting in long, drawn-out arguments between 
Manufacturing and Engineering. The only answer I know is for 
Management and Manufacturing to have confidence in Engineering, 
but at the same time people in Engineering must do enough work to 
be sure they are right. 
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INVOLUTE ANGLES GEAR 
Figure 7. Pinion & Gear Tooth Modification 


Since this investigation was so fruitful and can be applied to 
other situations, it is believed worthwhile to discuss the items just 
mentioned in detail. 


PROFILE 


A layout was made of the gears shown in Figure 4. The surface 
pressure was determined by the Hertz method, and the sliding 
velocity of the gear tooth surfaces calculated. The product of these 
two values was taken for each point of contact, under the assumption 
that only one mesh of teeth carries the load through the action, as 
the proper division of load between two sets of teeth in mesh is 
difficult to obtain. The method of calculation is shown in Figure 5. 
It is somewhat similar to that used by Almen in calculating PV 
values of rear-axle gears. Since the product of the pressure and a 
frictional coefficient and sliding velocity should represent a meas- 
ure of the power loss or heat dissipation for each instant of time, it 
was felt the action should be balanced so this product of pressure 
and sliding velocity is the same at each end of the action, The 
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Figure 8. PV Values for Two Sets of Teeth in Action 


frictional coefficient was assumed constant for each given value 
of pressure. The products of pressure and sliding velocity were 
designated PV values to agree with previously-established nomen- 
clature. The actual PV value is admittedly zero at each end of the 
action, but by balancing the values at the start and end of action, the 
ends of the single-mesh action will also be substantially balanced. 
Figure 6 shows the PV values for one set of teeth in action. The 
original and revised length of action is indicated, The original con- 
tact ratio of 1.68 reduces to 1.35. This indicates that contact ratio 
by itself is not a good measure of gear capacity. These figures and 
the test experience resulted in increasing the amount of metal re- 
moved from the form of the gear and the pinion teeth, as shown in 
Figure 7, This figure shows the pinion tooth modification above the 
zero line and the gear tooth modification below it. The action of the 
teeth moves from left to right. Clearance A (between respective 
contacting points at the root of the pinion and tip of the gear at the 
start of action) minus clearance B (between the contacting points of 
the teeth in the preceding mesh) represents the net clearance or 
“net relief? which is shown in the lower part of the chart. This 
relief prevents undue contact pressure which would result from the 
deflection of the teeth in the preceding mesh under load, The in- 
crease in the clearance or net relief is shown by the shaded area in 
the net relief chart. The changes shown in Figure 7 balanced the 
action of the gears so that scuffing tendencies should occur no more 
at the start of action than at the end. 
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Figure 9. Parallelisn f Gear Teeth 


Figure 8 shows an estimated curve of actual PV values con- 
sidering two sets of teeth in contact compared to the “single mesh” 
values for one set. 


MOUNTING 


The method given by Rolls-Royce for aligning the gears at the 
time Packard started production consisted of bluing the teeth in 
assembly to secure a heavy bearing at the rear of the teeth, by 
swinging the reduction gear cover in a rabbet concentric with the 
propeller shaft, and then doweling the cover in place. 

A careful check on approximately 200 engines indicated that 
scuffing appeared more frequently at the rear of the teeth when the 
bearing was heavy at the rear. Deflection tests of the gears in a 
four-square type of test fixture were carefully analyzed. They 
indicated that the axes of the gears assumed a new position equiva- 
lent to being transposed parallel to their original axes. This infor- 
mation resulted in doweling the reduction gear cover to the crank- 
case with the pinion bearing bores in line and parallel to the 
reduction gear bearing bores. Better bearing across the teeth was 
obtained from this procedure. 
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surface Roughness Measurement 


PARALLELISM 


When the gears were blued-in, a considerable variation in 
parallelism of the teeth was observed, clusters of five or six teeth 
being out of line. Very careful parallelism checks of every tooth in 
each of a dozen reduction gears were made after each step of the 
processing. A record of the checks on one gear is shown in 
Figure 9. 

It was found that most of the distortion occurred in carburizing, 
the gears going out-of-flat as well as out-of-round, A change in 
processing eliminated this trouble. It was also found that when the 
gears were too much out-of-parallel before grinding, it was difficult 
to remove all the error by grinding. The teeth of the reduction gear 
and the pinion were finally held parallel within 0.0003” in the width 
of the tooth without difficulty. Parallelism has been stressed, as it 
has been found much more difficult to maintain than proper index. 
Since little difficulty was experienced in maintaining index within 
0.0002” between the adjacent teeth of these gears, the effect of 
index errors on scuffing was not determined. 
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Figure 11. Reduction Gears After Model Test 
SURFACE FINISH 


In order to improve the grinding technique, the grinding wheels 
were dressed coarser with the dressing diamonds. Wheels dressed 
too fine tended to load up, particularly at one point, and give a 
pattern which appeared to be the result of a vertical bounce and 
horizontal oscillation, apparently induced by the exciting force of 
the loaded wheel. This pattern consisted of alternate high and low 
spots on the teeth, those on one side of the tooth space being dis- 
placed from those on the opposite side by the distance the table 
traveled during one-half wheel revolution. Coarser dressing elimi- 
nated this difficulty and also reduced magnaflux checks, particularly 
at the roots of the teeth. 

Observation of a large number of the gears after the acceptance 
run indicated the coarser finish showed less marking and “fussi- 
ness” than the fine finish, Comparison of the profilometer readings 
with profile recorder charts, as shown in Figure 10, led to specifying 
surfaces to be made in accordance with samples to obtain the char- 
acter of surface found desirable. Finishes, in agreement with the 
samples, were accepted between profilometer readings of 15 to 37 
root mean square. Under proper conditions, gears of this type of 
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Figure 12. Empirical Formula for Comparing Tooth Deflection 


of Comparable Gear Structure 


finish will go through a model test without wearing through the 
grinding marks. Figure 11 shows a set of gears after model-test. 

As a result of the improvement in gear tooth finish, the same 
method of finishing was applied to camshaft rocker levers. Revision 
of the finish specification from 10 maximum profilometer reading to 
new limits of 15 to 23 resulted in a distinct improvement after 
acceptance runs. 

This finish was applied to steel gears running with steel gears 
only. In one instance it was found necessary to make a crossed 
axis helical steel gear smoother to prevent it from cutting the 
mating bronze gear. 


LUBRICATION 


Improper lubrication at the start of the acceptance runs was 
suspected. A pump was installed on one dynamometer stand which 
was started when priming the engine, and was run until 10 minutes 
after starting under power. The record of scuffed gears of this 
stand for approximately 175 engines was compared to the general 
average of all the other stands, and a very definite improvement 
was noted. The acceptance test schedule was revised to call for 
priming lubrication of the reduction gears before starting the 
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engine, and continuing until 10 minutes after firing, to ensure suffi- 
cient oil to run in the gears properly. 


SUMMARY 


Careful attention to all the details listed resulted in reducing 
scuffing, from an objectionable percentage of rejections in the 
factory, to less than 0.5%. Service experience has shown that if the 
gears do not scuff during the acceptance run, they are not likely to 
scuff in service. 


NEW GEAR DESIGN 


Further anticipated increases in power indicated the desirability 
of more capacity in the reduction gears. A new design was studied 
with a pressure angle of 25 degrees and a pitch of 4.1. The length 
of action decided upon was based on permitting PV values of 
3,000,000 at each end, resulting in a contact ratio of 1.25. After 
establishing the usable action, the gear teeth were proportioned for 
strength. The modification was then established, using the previous 
deflection data of the production gears and applying the empirical 
formula shown in Figure 12. 

It will be noted that the teeth are assumed as deflecting about a 
point in the middle of the rim. Calculation of the deflection of a 
tooth as a non-uniform beam will account for only a small percent- 
age of the actual deflection present. The tooth actually bends as a 
body under load. The net relief shown below the modification chart 
equals A minus B as noted on the chart for start of action. A is 
taken where contact actually starts, and the net relief represents 
the clearance necessary to provide for the deflection of the teeth 
under load in the preceding mesh. 

Between the point of first contact and the outside radius of the 
gear, there is a ramp of 0.030” length on the gear tooth. Visual 
inspection of this ramp indicates if the contact is starting where 
desired. This determines if the proper amount of net relief has 
been provided. This also enables one to establish the amount of 
deflection in the preceding teeth by referring to the charts of the 
gears observed. The net relief at the end of action is established 
using the same procedure. 

The ramp at thetip is very important to assure that first contact 
does not extend to the tip of the gear tooth. The radii of curvature 
of both gear teeth should contact without loss of area to keep the 
unit pressure down and improve lubrication. 
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AXIS OF ROTATION 





STRESSCOAT INDICATIONS 


Figure 13. Reduction Gear Stresscoat Study 


The stresscoat checks, as shownin Figure 13, help to confirm 
the theory that the gear tooth deflects about a point in the center of 
the rim. 

Figure 14 shows the final modification charts for the gear and 
pinion as they appear on the respective engineering drawings. They 
have been matched together to enable the net relief to be determined. 
Separate charts were kept as a permanent record of each set of 
production gears. 

A closed circuit reduction gear fixture, shown diagrammatically 
in Figure 15, was designed and built to test the new gears at loads 
in excess of engine loads. Loading of the gears was accomplished 
by winding up a calibrated torsion shaft between the reduction gear 
pinions by the turnbuckles on one propeller shaft flange. The 
torsion shaft was enclosed and supported in a tube to provide lateral 
stiffness. On this fixture the test gears were lubricated and cooled 
by low-pressure jets of oil directed into the mesh, as in the engine. 
The dummy gears at the opposite end were lubricated by high- 
pressure oil jets on the teeth after mesh. The dummy gears were 
first lubricated by an e-p oil of SAE 90 viscosity to give them more 
surface capacity than the test gears. Later, engine oil was used, 
but in greater quantities than on the test gears. 

Test running of the new gears indicated their capacity to be 
appreciably greater than that of the original gears. The action 
started and ended where called for by the design, with no scuffing 
except when overloaded appreciably, with the result that no change 
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whatever was made in the modification established in the original 
design. A year and a half of production experience with no scuffing not ; 
indicated no change necessary. The modification formula has been gear 
used on several other gear-sets and found satisfactory, but is woul 
offered with hesitance since a large number of different gear load 
designs have not been checked. Care must be used in obtaining the coef 
original deflection data. The designs compared should have the the | 
same general structure. are 
The closed circuit test fixture gave an opportunity to determine 
the efficiency of the gears, which was approximately 99% per mesh. to a 
extr 
ticu 
PLANETARY GEARING the 
tion 
Planetary gearing, used in overdrives inthe thirties, has come mat 
into prominence in automatic transmissions, particularly since the 
World War II. These gear sets will carry large powers for their exc 
sizes because the load is divided between several meshes of gears. pitt 
Alignment of the teeth is very favorable, with regard to deflections wit! 
of mountings, and the product of the loads and the pitch line ve- | __ witl 
locities (or more properly the base circle velocities) is ordinarily whe 
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WIND UP MECHANISM 
TO APPLY TORQUE TO 
THE GEAR SETS 


SHAFT CONNECTING GEARS TEST GEAR 
ENGINE 
DUMMY GEAR enaamernan 
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TO DRIVING TORSION 
DYNAMOMETER SHAFT 
CONNECTING 
DUMMY PINION PINIONS 
Figure 15. Diagram of Reduction Gear Test Fixture 


not as large for a given reduction as in multiple reduction layshaft 
gearing in series. Strictly speaking, the power loss inthe gearing 
would not be measured directly by base radius velocities-times- 
load-times-a factor, but the load on the teeth-times-a “suitable” 
coefficient-times-the sliding velocity on the teeth summed up over 
the line of action. Determination of load on the teeth where two sets 
are in mesh is, of course, not easy. 

Reduction of power loss is as important in keeping the teeth cool 
to avoid distress as to save fuel. Proper lubrication of teeth is 
extremely important. This applied also to needle bearings, par- 
ticularly those supporting the pinions. I recall one case in which 
the gear-set, with marginal lubrication as furnished in the produc- 
tion unit, permitted the pinion needle bearings to weld in approxi- 
mately 30 minutes under full engine torque. Adequate lubrication to 
the bearings enabled the gears to run 50 hours without distress 
except pitting of the pins on which the needle bearings ran, This 
pitting started rather early, but we continued to run the gear-set 
without any ill effects until considerable pitting had occurred; but 
with the set still entirely operable. That raises the question as to 
when we should classify a unit as “failed” due to pitting. 
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Modification of the involute is desirable on planetary gearing 
but from the writer’s experience, appears to be only enough to} 
prevent interference on start of action; as the deflections appear ty 
be less than those on gear trains like aircraft reduction gears 4; | 
described. In any case, one can look at the teeth and see wher | 
contact starts. That is the final answer. 

Helical contact ratios of the order of 1.3 to 1.4 have resulted jy | 
satisfactory quietness in gear sets with which the writer is familiar, 
These ratios can be obtained due to the finer teeth, with helix angles 
of the order of 18 to 20 degrees, without excessive gear width. This 
makes axial thrust and gear tooth lead alignment problems easier, 
Very slight crowing of teeth and very slight cutting of pinion teeth 
“off lead” has been attempted at times to maintain proper bearing 
on the teeth. 


CONC LUSION 
To make highly-loaded gearing work successfully, we should: 


1, Make mountings as stiff and as accurate as possible to 
secure accurate alignment of tooth faces. If there is flexibility, it 
should be of a nature that causes the teeth to align and not misalign 
themselves. 

2. Make gearing accurate for lead or parallelism. To be suc- 
cessful, gearing must be held more accurately than the production 
people are willing to see shown on the drawings. 


3. Modify the teeth properly to suit the conditions to which they | ge 
are subject. An analysis will show you can run on only about sc sé 
much of the profile due to excessive sliding if you go too far under a 
high loading. nm 

4. Secure proper surface finish to suit the conditions present. m 

5. Provide proper lubrication, particularly at start of operation, fa 

6. Secure the best balance between surface durability and beam tc 
strength of the teeth in the original design. We make teeth fine t vi 
obtain surface durability; we make them coarse to obtain beam a 
strength, Breakage of teeth is, of course, a greater hazard than | d 


scuffing or scoring them, but gearing has been over-designed in this 
regard at the expense of surface durability. 

7. PV values as described in this paper or PVT values as de- | 
veloped by Almen and Straub are the best methods the writer knows 
for determining surface capacity of gear teeth. PVT values were : 
developed to overcome one defect of PV values which showed up | 
under large variation of pitch line speeds and sliding. We used PV 
values because of our experience background. 
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Figure 16. Major Gears-Rolls-Royce Aircraft Engine 


By observing the rules just listed, we were able to make all the 
gears in the previously-mentioned aircraft engine operate very 
satisfactorily under high loads or high speeds (Figure 16 shows a 
set of the major gears). These gears were not quiet, as it was not 
necessary to make gears run quietly under high loads. I feel it is 
much more difficult to make them run quietly than just run satis- 
factorily. This is the problem we have in the units in the automobile 
today. Making them run is not enough — they must also be quiet and 
vibration-free. This is often a difficult problem; however, it is 
another subject, and a big one which cannot be included in this 
discussion. 
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